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A great variety of experiments, like heat release measurements, acoustic measurements, and trans- 
port measurements on mesoscopic samples have proved that two level systems (TLSs) have a crucial 
role in the low temperature thermal and electric properties of disordered systems. This paper is 
aimed at reviewing the role of slow TLSs in point contacts. First the theory of point contacts is 
summarized, concentrating on the discussion of different point contact models, and on the different 
regimes of electron flow in the contact, mainly focusing on the ballistic and diffusive limit. The 
Boltzmann equation is solved in both regimes, and the position dependence of the electrical po- 
tential is determined. Then the scattering processes in point contacts are investigated, particularly 
concentrating on the scattering on slow TLSs. If the the electron assisted transitions between the 
two states are negligible the electron-two level system interaction can be treated with a simplified 
Hamiltonian. The scattering on such slow TLSs causes nonlinearity in the current-voltage charac- 
teristics of the point contact, which can be determined using Fermi's golden role. These calculations 
are presented showing both the contribution of elastic and inelastic scattering, and including the 
dependence on the position of the TLS, and on the effect of high frequency irradiation. These 
results are used to discuss the differences between these slow TLSs and the fast centers which may 
be described by the two channel Kondo model. The available experimental results are analyzed, 
distinguishing between the effects due to the different types of TLSs. 
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I. INTRODUCTION 

In the last decades point contact spectros copy has been proved to be a very pow erful method to study electron- 
phonon interaction in metals (for reviews see lYansonl986alJansenl980l and lDuiti989>l . In the recent years it has been 
realized that in mesoscopic metallic systems the electron scattering by some dynamical defects may play a major role 
in unders tand ing the energ y relaxation and dephasing processes of the electrons at low temperature (for reviews see 
lPecs200(]| and IUisaghv200ir) . Such defects can also be studied by point contact spectroscopy provided that they are 
situated in the contact region. These defects can either be magnetic impurities or structural defects where an atom 
or a group of atoms have a metastable, almost double degenerate groundstate forming a two- level system (TLS). The 
TLS spectroscopy is very similar to the phonon spectroscopy, but as the TLSs cannot propagate like phonons, they 
are not in thermal equilibrium with the bath. The TLS can be formed as defects in the bulk region of crystalline 
materials, or at dislocations, at surfaces due to impurities and also in amorphous materials. In most of the cases the 
internal structure of the TLSs are unknown and they can behave very differently. 

The dynamics of magnetic impurities shows the Kondo effect^^SS^SiSS where by lowering the temperature a Kondo 
resonance is developed at the range of the characteristic Kondo temperature, Tx and approaching zero temperature 
a magnetic singlet is formed by the impurity spin and the electrons. In that case near zero temperature the magnetic 
singlet acts just like a static impurity and the dynamical processes are frozen out. 

The subject of slow defects can be described by the conventional theory of TLSs and its application to the point 
contact is the present review. There are, however, many experimental facts which cannot be described in that way, 
and show resemblance to the magnetic Kondo anomalies even if their magnetic origin is very unlikely. We make an 
effort to point at those experiments by comparing the observed behavior with the predictions of the slow TLS model. 
Several of these are believed to be due to dynamical defects, different from slow TLS. It must be emphasized that there 
is no theory which is generally accepted and can be applied. Earlier it has been proposed that they are fast TLSs, 
which can be described by the so called two-channel Kondo modeli^^'^^'^''^^^'^'^i^''^'^'^''^^^'^'^i^''^'^^''^^'^^'= That proposal has 
recently been criticized and debated.^^"'^'^''^°°^'^A''''"^''^°°^'^i^i"'^R^Y200iiBorda2003 r^Yiese scatterers now appear rather as 

"phantoms" and the present review does not cover that issue, we, however, present the earlier interpretation just to 
help in reading the earlier experimental publications. 

(i) slow TLS . In that case an atom or a group of atoms can have a transition between the two states of lowest 
energies by direct tunneling or at higher temperature by thermal excitation over the barrier. In the former 
case those systems are also called as tunneling centers (for review see2i2£Jii22i) . The transition rate can be on 
the scale of seconds and even on longer scales. In that case e.g. the specific heat may depend on the speed of 
the measurement as to reach thermal equilibrium takes longer time. The resistivity of such mesoscopic system 
may jump in time between two or several values on the above time scale, as in small mesoscopic systems the 
resistivity depends on the position of a single atom even if the two possible atomic positions are much smaller 
than an atomic distance. The study of point contacts containing TLSs in the contact region turned out to be 
very powerful as the TLS may result in zero bias anomalies shown by the I{V) characteristics. The electron 
passing through the point contact can interact dynamically with the TLS or it can be elastically scattered, which 
scattering depends on the actual position of the atom forming the TLS. The inelastic scattering may result in a 
back-scattering of the electron in which case the electron does not proceed from one electrode to the other when 
the applied voltage is larger than the splitting of the TLS and the process shows up as a resistivity minimum at 
zero bias. In the case of elastic scattering the transmission rate may depend on the position of a single atom and 
the anomalies can have either sign. The characteristic relaxation times from the TLS causing slow telegraph 
fluctuation of resistance to those being responsible for zero bias anomalies in the I{V) curve varies over several 
orders of magnitude, but for all these so-called slow TLSs the average transition time is much longer than the 
electron-TLS interaction time. 

(ii) fast centers. The fast TLSs are strongly debated whether they can show Kondo 
anomalyj ^^''''""'^W^^'Vi'"^S^Y^W^i^°'''^''^W? Here we present the previous ideas just to make contact with the 
extended literature, and the present status of the problem will be discussed at the end of the paper. If a 
conduction electron experiences several switching occasions one after the other we speak about fast TLS. These 
fast TLSs may show Kondo anomaly as the TLS's two states can be described by a quasi-spin replacing the real 
impurity spin of the magnetic Kondo problem and the scattered electron can have different spherical momenta 
e.g. s- and p-states playing similar role to the real spin of the conduction electron. The spin-flip process of the 
magnetic Kondo problem is replaced in the latter case by electron assisted transition between the two states of 
the TLS. That trasnition was is originally believed to be electron assisted tunneling, which is unlikely according 
to the present theoriesj^iSiSSSSllflia^iSiiiaSfilii^ however an electron assisted transition without tunneling is still 
possible i^2£^S222^ In point contacts containing magnetic impurities the Kondo resonance scattering contributes 
to the back-scattering rate, thus the Kondo impurities result in zero bias peak in the dynamical resistivity. 
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which has a width characterized by the Kondo temperature, Tk- The sign of the peak is just opposite to those 
caused by slow TLSs. I.K. Yanson and his co-^orVeva^'''^''°'^^'^^^''^^y''^^^'^'^^°P'^''^y'''^^^^^ 

and also other groupsiisagaiiSSSiSaiEiiiSS studied the zero bias anomalies due to intentionally placed magnetic 
impurities. Considering the Kondo effect there is a significant difference between the magnetic case and the 
TLS. In the latter case the conduction electrons, additionally to their angular momenta occurring in the 
Hamiltonian and the coupling, have a further degree of freedom, the real spin, which makes the electronic sea 
double degenerate although it docs not appear in the coupling constants. Referring to that extra degeneracy 
the phenomenon is called two-channel Kondo problem (2-CK) in contrast to the magnetic one-channel case 
(1-CK) (for earlier review see ICoxl999ll . That extra degeneracy prevents the system to have a singlet ground 
state at low temperature, it has finite entropy and, therefore, the Fermi liquid behavior does not show up as far 
as a low energy cutoff (e.g. the energy difference between the two lowest states of the TLS, known as splitting) 
is reached as the temperature is lowered. 

The condition for the formation of a 2-CK ground state is that the energy splitting is negligible com- 
pared to the Kondo temperature, which is generally not expected. The zero bias anomalies have been 
believed to be powerful t ools to make difference between 1-CK and 2-CK problems (for a review see 
lvonDelftl998lvonDelftl999^ . For example annealing or electromigration could result in disappearance or modifi- 
cation of the zero bias anomalies in the case of structural defects while such effects are not expected for magnetic 
impurities. Similar issues also appear in the electron dephasing time which has been intensively de- 
bated recently^^'T.i^^AW"'^ii999'A'""°'^200ia.Aicincr200ib.Uisap;hY200i,imrYi999_ Another version of the two-channel Kondo 

effect was proposed to explain the zero bias anomalies in URu2Si2 point contactsS2^£i6Si22i where two degen- 
erate localized electron orbitals of the U atom play the role of the TLS which interact with the conduction 
electrons, c°"i987,Coxi999 

The issue of the 2-CK problem was responsible for the great interest in the zero bias anomalies due to TLS and, 
therefore, it has a great importance to distinguish between anomalies due to slow TLSs and fast scatterers. 

Even if a considerable amount of the discussions are carried out considering fast TLSs there is no doubt that in 
many cases the slow ones are dominant e.g. where the zero bias anomaly in the point contact spectrum has a positive 
g^g^ Akimenkoi993,Keijsersi995^Keijsersi998 forking ou the issucs listcd above we learned that the literature on the theory 
of slow TLSs in point contacts is very much scattered over many papers and journals, the predictions should be 
collected from a large number of papers, especially, published in Russian and Ukrainian physical journals some of 
them not available at smaller libraries. That experience encouraged us to collect and at several occasions to complete 
those results. The present review is written in a self-contained manner, thus no further reading is necessary to follow 
the theory. 

We want to emphasize that the authors did not consider as a task to give the historically proper complete list of 
references but they picked up those which are the most appropriate ones to get further information. 

II. TWO LEVEL SYSTEMS IN SOLIDS 

A. The model of two level systems 

The model of two level tunneling systems was at first introduced to explain the unexpected low-temperature behavior 
of specific heat and thermal conductivity in disordered solids which cannot be described by phonon excitationsi^£ii££i22i 
The model based on two level systems, which was able to explain the linear temperature dependence of specific 
heat of amorphous dielectrics at low temperatures was suggested by Anderson, Halperin and Varma^S^SSSSiSS and 
Phillips£iiiiiifiSi2S This so-called tunneling model assumes that in an amorphous state some atoms or groups of atoms 
may switch between two, energetically nearly equivalent configurations. This situation is modeled as a double well 
potential with two stable states differing by energy A, and a potential barrier between them with a tunneling rate Aq 
(see Fig.^. The tunneling probability can be approximated in terms of the Gamow parameter (A): Ag ~ e~^, where 
A is defined by the characteristic height, V and width, W of the potential barrier, and the mass, M of the tunneling 
particle: 



Furthermore it is assumed that the TLSs in a solid are uniformly distributed in terms of the parameters A and A, 
that is the density of states P for the first excited states of the TLSs is constant: 

P(A, A)dAdA = PodAdA. (2) 




^ ► 

d 

FIG. 1: A two-level system is modelled as a double well potential with an energy difference A between the two positions, and 
a tunneling probability Ao for crossing the barrier between the two metastable states. W and d denote the width of the barrier 
and the distance between the minima, respectively. 



Experimental studies do indeed support this assumption."""''""S'^''i?7?,9°idinfii^7?,Biacki?^i 

The two positions of the system are described by a quasispin, and the Hamiltonian of the TLS is expressed in terms 
of the quasispin's Pauh-operators: 

Htls = ^ E i^b+a^^ph + ^ob+a^pbp) , (3) 

a, (3 

where the indexes a and f3 correspond to the two states of the TLS (|1) and |2)), and and ba are the creation and 
annihilation operators for the state a. This Hamiltonian can be diagonalized by writing the eigenstates as: 

\E+) = fi\l)+iy\2) 

\E^) = H1)-m|2). (4) 

The eigenenergies are: 



E^^±l^A^ + Al (5) 



and thus the energy splitting between the two eigenstates is: 



E = ^A^ + Al. (6) 

The coefficients /x and v can be expressed as: 

E + A 



^{E + Ay+Al 

Aq 

^{eTWT^o' 



(7) 



For a highly asymmetric TLS, where A 3> Aq the energy eigenstates are the two positions of the TLS, that is = 1, 
and v = 0,1. In the opposite case, where the tunneling energy is much larger than the splitting (Aq S> A) the energy 
eigenstates are delocalized between the two positions, and ii = v = 



B. Scattering on two level systems 



In a disordered material two level systems are not isolated objects, but they interact with their neighborhood. For 
instance a strain induced by an elastic wave may change the parameters of the two level systems. It is assumed that 
only the asymmetry parameter A is influenced by an elastic wave, thus the interaction between two level systems and 
phonons can be treated as a perturbation in A: 

^A = -/26u^k, (8) 



where 7 is the coupling constant, and Suik is the appropriate local elastic strain tensor. 
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In metallic samples the coupling to conduction electrons has to be taken into account as well. The general form of 
the Hamiltonian describing the interaction between the electrons and the TLS isiZawadowskii980a,Zawadowskii980b 

pi,P2,Q!./3,(T 



where a+ and a are the electron creation and annihilation operators, a is the spin index of the electron and V is the 
interaction matrix element, which can be decomposed in terms of Pauli operatorst2Sfi2^S21§JSii22iik 



PlP2 



E 

=x,y,z 



Pa 



(10) 



The last term stands for potential scattering on the average positioned TLS, the term with describes the difference 
between scattering amplitudes for the two positions of the TLS, while cr^ and correspond to the electron assisted 
transitions between the two states of the TLS. 



1. Slow two level systems 



The interaction Hamiltonian is further simplified if the electron-TLS interaction time is much shorter than the 
transition time of the TLS, thus the TLS stays in a given position during the electron scattering. For these so-called 
slow two level systems the electron assisted transitions described by and can be neglected, and the Hamiltonian 
is: 

i?ci-TLS = ^PiP2«P2<Tapia E ^p"iP2«^2<Tapi'T (11) 

the interaction matrix elements between the delocahzed energy eigenstates are: 

{E^\H,i^TLs\E+) = [2/^^^P2pJ «k^«Pi- (12) 

Pl,P2,0- 

which gives the amplitude of the electron's inelastic scattering on the TLS, which is associated with an energy change, 
E. The matrix elements corresponding to the elastic processes are written as: 

{E+\H,i^TLs\E+) = Y [Cpi + (a^' - '^')^P2pJ «P2.«Pi- 

pi,P2,<T 

{E.\H,i.TLs\E^) = E [Cpi - - ^XV] «P2-«Pi- (13) 

pi,P2,cr 

In both cases the energy of the scattered electron is conserved, however the scattering cross section is different if the 
TLS is in the or \E^) state. 

By simple argumentation the interaction strength, can be estimated as: ~ (^p^^-^y Biax^kidrB^Yiada.iWSSa ^ 
the difference between the electron scattering amplitudes in the two positions must be proportional to the distance 
of those positions on the scale of the Fermi wavelength, and the strength of the scattering, V. V'^ is very model 
dependent, e.g. the tunneling atom can be different from the most of the atoms. 



2. Fast two level systems 



The electron assisted transition processes are only important in the case of fastly relaxing two level fluctuators. In 
this case the Hamiltonian given by Eq. ^ can be scaled to the two channel Kondo problem (2CK) with the Hamilton 
operator: 

Hint ~ (TTLSV ^ Gp^^^SalSap^fSa, (14) 

pi,P2,Ct,/3,(T 

where the the momentum depe ndence of On n is decompose d into appropriately choosen spherical waves with indices 
a, (3, and p = |p| (For details see lVladarl983alVladarl983b(l . S is a pseudo-spin corresponding to the spherical indices 
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of the electrons, whereas <ttls is a pseudo-spin corresponding to the states of the TLS, and V is the interaction 
strength. In this model the electron spin is an extra degree of freedom compared to the single channel magnetic impu- 
rity Kondo problem, and this quantum number is conserved during the interaction and the coupling is independent of 
that. In contrast to the single-channel Kondo problem, which has a non-degenerated, S = ground state because of 
the spin compensation cloud formed by the conduction electrons, the ground state of the two channel Kondo problem 
has a non-zero entropy. In the two channel Kondo model the conduction electrons overscreen the impurity spin by 
forming a non-trivial spin state. The low temperature regime of the two channel Kondo model was found to have 
non- Fermi- liquid behavior. 

This paper is devoted to review the literature of the slow TLSs and point contacts. The Kondo regime of the 
TLSs is out of present consideration, but one can find an up-to-date review of this field in the work of Cox and 
Zawadowski£2£i22S 



3. Scattering rate beyond the Born approximation 

In the previous part of this Section the scattering rates due to TLSs are calculated in Born approximation. Con- 
sidering only the screening interaction V and ignoring the assisted tunneling rates Kondo^^sa^aiSSi has 
found power law correction to the tunneling rate which is due to the building up of a screening cloud or Friedel 
oscillation around the changed position of the tunneling atom. That is closely related to th e Anderson's orthog- 
onality catastropho-^'^'^'^''"™^?^7iK°ndoiii)7^,Kondoit)^^ ^^^^ ^-^e X-ray absorption problem (see e.g. lNozieresl96fl() . That 
phenomenon is associated with creating a large number of electron-hole pairs as the position of the tunneling 
atom is changed. The renormalized dimensionless coupling constant (usually denoted by K) can be expressed in 
terms of the phase shift S describing the scattering of the electrons by the atom in the s-wave channel (see e.g. 
lYamadal985lKondol988j) and the separation distance, d between the two positions of the TLS 



K 



1 Vl — tan 5 

-arctan^^_^^2^^^2^^i/2 



(15) 



where a = jo{kpd) with the spherical Bessel function, jo and K < 1/2 for only s-wave scattering. In the weak coupling 
limit (5^1 and fcp d <g; 1 that expression is proportional to {gaV'^)'^. In many other publications the notation a is 
used for K (see e.g. ICostil996|) . 

The behavi or of a TLS co upl ed to electronic heat bath (Ohmic region) by the coupling K has been studied in great 
detail fsee e.g. lLeggettl987l and lWeissl9991 In general there are two different regions regarding the coupling strength 
K. Here the asymmetry parameter A is disregarded and T = is taken. 

(i) < K < ^ damped coherent oscillation. At K = the atom has a periodic motion between the two positions. 
As K is increased the oscillation is more damped and the renormalized tunneling frequency Aren is given by 



Aren — Aq 



Ao 



(ii) 



where lUc is the high frequency cutoff, which has first believed to be the electronic bandwidth (see e.g. lCoxl999t 
but according to the adiabatic renormalization it is rather the energy scale of the lower excitations in the 
potential well where several excited levels are considered i^SSSiiSSSil^Sa^aiSS 

^ < K < 1 incoherent relaxation. The motion is incoherent and the atom is spending longer times in a position 



before tunneling as K is increased. 



Many experimental facts indicate (see next subsection) that the TLSs are in the first region {K < 1/2) and even can 
be in the weak coupling limit {K <C 1). 

The main issue is how the quasielastic and inelastic scattering rates are modified by the coupling K. That problem 
was investigated by H. Grabert, S. Linkwitz, S. Dattagupta, U. WeisaSSi^££ii22& considering the neutron scattering 
by TLS. Those results are very instructive for the present case, there are, however, essential differences: 

(i) Korringa type of relaxations are very important in the relaxation of the TLS coupled to the ohmic electronic 
heat bath. If energy uj is transferred to the TLS the spectrum of the TLS is broadened by a relaxation rate 
proportional to lu. In the case of point contacts that energy may arise from the energy of the electron passing 
through the point contact. At finite bias, V the system is, however, out of equilibrium and in the stationary 
state of the TLS the averaged energy of the TLS is proportional to V also. Thus, the nonequlibrium situation 
should be treated by using e.g. the Keldysh technique i^Si^SSiiiSSi 
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(ii) The energy transfer is not determined by the applied voltage, in contrast to the neutron scattering experiment 
where the energy of the incoming and outgoing neutrons are simultaneously measured. In the present case by 
small change of the voltage the energy of the extra electrons passing through the point contact are given (~ V) 
but the energy of the scattered electrons is not determined, thus the energy of the outgoing electron shows some 
distribution. That results in an additional smearing in the spectrum. 

In the following we give a brief summary of the neutron scattering case£Sik£££i22& Very accurate results are ob- 
tained for low temperature and energies except the range u <^ Are„. That is certainly the interesting range con- 
sidering especially for the long range tails of the characteristics which are crucial for the interpretations of the 
experimentsi^^SlSgiSiSSfei^SS"'^!^^* The calculation was carried out by expanding in the tunneling events and summing 
up. 

More comprehensive studies were performed by using numerical renormalization group 
energy range of energy and even including the asymmetry, A and later even for nonequilibrium£2§iiiS2i. 
The previous resultsSs^££ii22& have a very good fit, by the formula 

b[UJ 



exp(fcr) - 1 (w - ujo)^ + huj^ 



where a, B, and luq are fitted parameters. That clearly shows a double peak structure with energies iwo. The lines 
are broadened by a Korringa type energy relaxation rate proportional to a;. In the nearly Lorentzian lineshape no 
direct role of coupling K was discovered. 

Those results are strongly indicating that in weak coupling region that would be very hard to discover any sign 
of anomalous relaxation rate in the tunneling characteristics. That situation is further complicated by the averaging 
over the distribution of the energy of the scattered electrons and the out of equilibrium situation. We find, that the 
experimental studies of the tunneling characteristics are not appropriate tool to look for the interaction effects, except 
that Aq is renormalized. That situation is very different from the cases of one and two channel Kondo impurities, 
where the characteristics are determined by the Kondo resonances. 



C. Experimental investigation of two level systems 

There are several experimental methods to investigate the kinetics of TLSs. The first two methods to be discussed 
are studying the TLSs by their interaction with phonons. In heat release measurements the time dependence of heat 
release or specific heat is measured: after a sudden cool down of the sample due to the energy relaxation of the 
TLSs it takes a long period of time to reach thermal equilibrium. The second group of measurements investigates the 
propagation of sound in amorphous materials measuring sound attenuation and sound velocityiSliSisiiiiS££i22& 

Heat release measurements basically focus on amorphous dielectrics and superconductors, as in metals the con- 
tribution of conduction electrons to the heat capacity usually dominates that of TLSs, thus the estimation of Pq is 
more difficult. According to the heat release measurements of Kolac et al>J^Si2£i22& in amorphous metal Fe8oBi4Si6 
the density of the TLSs is about Pq ~ 1-5 x IQ'^^J-^m-'^. 

Though in acoustic measurements basically the phonons are investigated, if the electron - TLS coupling is strong, 
the results may significantly differ from that on amorphous dielectrics £si^iiiSi22& A simple description for ultrasound 
absorption by two level systems is as follows. The occupation numbers for the upper and the lower states of the TLS 
are n+ and n_ = 1 — n^, respectively. The probability for ultrasound absorption is aPri-, where P is the ultrasound 
intensity, and a is the absorption coefficient. The probability of relaxation is ti^/t, where r is the relaxation time of 
the TLS. In balance these two quantities are equal, thus one obtains: 



1 1 



T aP 



(16) 



The absorbed energy is: 



En 



hioaPn^ = hiu- 



aP^ 



aP 



(17) 



It is easy to see, that at large enough intensities (P — s- oo), the energy absorption saturates to a value inversely 
proportional to the relaxation time (Eats — > const. /t). In insulators and semiconductors the available maximal 
ultrasound power is enough to drive the absorption into saturation at T = 1 — 4K, but in metals such saturation is 
only obtained at temperatures as low as T = lOmK. It implies that the relaxation times in metallic samples are much 
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shorter than in disordered insulators. The interpretation for this experimental observation is, that the electron-TLS 
scattering processes are dominating in the relaxation of TLSs. Assuming a standard Korringa-like electron - TLS 
relaxation, where an electron-hole pair is created, at low enough temperatures one obtains: 

- - KPokT, (18) 

T 

that is, the ultrasound absorption measurements give the opportunity to make an estimation for the electron - TLS 
couphng parameter (K), and the density of the TLSs {Pq)- In fact the acoustic properties of disordered metals do not 
show universal behavior as in the case of insulating glasses. There are several theoretical approaches (the mentioned 
standard Korringa-hke relaxation,^''''^'*^^*^ the strong couphng theory of Vladar and Zawadowski,^'^^^''^^*'^^-^''*'^^''^"^^- 
the electron-polaron effect considered by Kagan and Pi-okof'e\iJ^SSaiii22Sil^2ii^2i2Si) that may explain the aspects of 
the relaxation of TLSs due to conduction electrons 

Bezup,iYi20oo gyjj^ taking into account the simple electron-TLS 
coupling given by Eq. jnj, the acoustic measurements can give estimations for the density of two level systems 
(Pq) or the electron-TLS coupling constant {K). In case of superconductors the Korringa relaxation is suppressed 
by the presence of the superconducting gap, but the gap disappears if a high enough magnetic field is applied on 
the sample, switching on the TLS-electron interaction£2Efi£SiSiHii22^ According to the measurements of Esquinazi 
et ali£saiii£2Sii22£ in the normal and superconducting state of PdsoZrro compared with the theory of Kagan and 
ProkofeAiJ^SSaaiSS^ the coupling parameter is approximately K ^ 0.4. The acoustic measurements of Coppersmith 
and Goldin g^°PP'"''^™'*^^^^^^ on the normal conducting amorphous metal Pdo.775Sio.i65Cuo.06 estimated the couphng 
constant as if « 0.2. For the density of the two level systems a lot of results are available in the literature, and it 
can be stated generally, that it is the same within a factor of less than 10 over a large variety of disordered material, 
being metallic glasses or dielectric amorphous material. (For TLSs with energy splitting less than 1 Kelvin it is 
approximately 1-lOppm.) A detailed review of heat release and acoustic measurements in disordered material can be 
found in .Nittkel9 98 Esquinazil99_8, and for earlier data see Vladarl983a 

Electron-TLS interaction can be studied directly by measuring electric transport in disordered metals. Point contact 
spectroscopy offers the possibility to investigate the properties of even a single two level system centered in the vicinity 
of the contact. The most spectacular sign of two level systems in point contacts is the so-called telegraph noise: the 
resistance of the contact is fluctuating between two (or more) discrete values on the timescale of seconds. One can 
estimate the average lifetimes in the excited (re) and ground state (r^) by recording hundreds of transitions, and 
fitting the resulting histograms of lifetimes {Pe{t) and Pg{t)) with exponential decay functions {Pe{t) ~ e^*/'^= and 
Pg{t) ~ e~*/'^9, respectively). The two life times are related to each other by the detailed balance, Te/rg = e-^/^bT ^ 
where A is the asymmetry parameter. If the TLS jumps to the other state, the electron screening cloud also needs 
to rearrange, which makes the jumps of the TLS slower. (The building up of the electronic screening cloud is related 
to a process which is similar to the X-ray absorption in metallic systems£i25i£££Si2S2il^Sa^2i2S, and which is also called 
as electron-polaron effect ii^SSaiii222il^asaii22S) From the theoretical point of view this slow-down of the TLS motion is 
treated as a renormalization of the tunneling amplitude ii^2adsi224 

/An\ ^ 

A,e„=Ao[^j , (19) 

where lOc is the bath cutoff frequency for which originally the electron bandwidth had been takenji^sadsi22i but later it 
was proposed that it is replaced by a typical energy of the next higher excitation in the potential wellii^SSgiiii222il^2adSi22& 
Golding et al£syiasi222( studied the two level fluctuation in mesoscopic disordered Bi samples. They argued that their 
samples were in the strong coupling limit, where the coupling of the TLS to conduction electrons is very strong 
compared to the tunneling matrix element, thus the latter can be treated as a small perturbation. In this limit the 
tunneling of the TLS is incoherent, because of the dephasing due to the fast oscillations of the electron bath. Here a 
theoretical scaling function can be set up: T^~^^jf = fK.Ao{f^BT/ A), which does indeed agree with the experimental 
result. (The energy splitting of the the TLS can be tuned by changing the electron density, which can be reached by 
applying weak magnetic fieldi^iSiiag£SBiii22i) In the limit kgT/ A ^ 1 the scaling function is constant, i.e. ■^f ~ j^^K-i^ 
that is the coupling parameter can be determined from a single fit, giving K « 0.24 for the particular TLS measured. 

The effect of fastly relaxing TLS in point contacts cannot be resolved as a telegraph fluctuation of the resistance, 
it causes an anomalous behavior in the voltage dependence of the differential resistance around zero bias, a so-called 
zero bias anomaly (ZBA), which will be a subject of detailed discussion in this review. 



III. POINT CONTACTS 



Metallic systems, where two macroscopic electrodes are connected via a contact with small cross section are called 
■point contacts" (PC) in general, regardless of the actual size of the contact area. It can be generally stated, that the 
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FIG. 2: The four main techniques for establishing point contacts, (a) Mechanically controllable break junction, (b) nanocon- 
striction, (c) spear-anvil geometry, (d) touching two wires. 



resistance of a point contact is mostly determined by the narrow neighborhood of the junction; therefore, a PC acts 
as a "microscope" magnifying all kinds of phenomena occurring in the small contact region. 

Several methods have been worked out to produce extremely small contacts between two conducti ng leads. Figure^ 
presents the four most important experimental techniques. The first one (Fig. Ill; lMullerl992(l is referred to as 
mechanically controllable break junction technique (MCBJ). Here the sample (practically a piece of metallic wire) is 
fixed on the top of a flexible beam, and a small notch is established between the anchoring points. The contact is 
created in situ at low temperature by breaking the sample with bending the beam, thus one obtains clear and adjustable 
junction on atomic length scale. The second method (Fig. [21d; RalDhl995j) uses nanolithography to establish a small 
hole in a silicon nitride insulating membrane by etching. If the etching is stopped just when the hole breaks through, 
the diameter at the bottom edge remains extremely small (~ 3nm) , well below the usual resolution of lithography (« 
40nm). In the next step metal is evaporated on both sides, creating a high-quality point contact device. This method 
can provide extremely stable and clear contacts on atomic length scale, however the contact diameter cannot be varied 
during the measurement. The third drawing (Fig. 12;) shows a similar arrangement to a scanning-tunneling microscope, 
a so-called spear-anvil geometry: a vertically movable, sharply tapered needle is pressed onto a flat surface. Finally, 
Fig. 1211 shows a simple technique, where the edges of two wires are brought into contact. In arrangements (a), (c) 
and (d) usually a differential screw mechanism is used to adjust the contact supplemented with a piezo-crystal for 
fine tuning. 

The first application of point contacts was carried out by Igor YansonSiS22iii224 to investigate electron-phonon 
scattering in nano junctions. He found, that the point contact spectrum, obtained as the second derivative of current 
with respect to the bias voltage {(Pl/dV'^) contains structure due to the electron-phonon interaction described by the 
Ehashberg function a2(w)F(w).l^li2iJssi£iii22^ This simple method for measuring electron-phonon interaction spectrum 
became a popular application of PC spectroscopy, however, it can be used to probe other electron scattering processes 
as well, like electron-TLS interaction, which is the central topic of this paper. 

Theoretically, point-contacts are considered as two bulk electrodes connected through a narrow constriction. The 
simplest, and most commonly used PC model is presented in Fig. EJi. This so-called opening-like point contact is 
an orifice with diameter d in an infinite isolating plane between the two electrodes. Another extreme limit is the 
channel-like PC: a long, narrow neck between the bulk regions with the length being much larger than the diameter, 
d (Fig.|2b)- The crossover between the both cases can be obtained by considering the point-contacts as rotational 
hyperboloids with different opening angles (Fig. O;). In most of the cases the shape of the PC does not influence 
the character of physical processes in the constriction radically, and the main parameter is the ratio of the contact 
diameter (d) and other characteristic length scales in the system. Three fundamental length scales are the mean free 
paths connected to different scattering processes (Z); the Fermi wavelength of electrons (Af); and the atomic diameter 

(dat). 

If the Fermi wavelength or the atomic diameter becomes comparable to the contact size, we speak about quantum 
point contacts and atomic- sized point contacts, respectively. These systems are reviewed in detail in .Agrait20 03. In 
this paper we consider contacts, that are neither atomic nor quantum, but that are small enough compared to certain 
mean free paths, and the calculations are basically performed for an opening type geometry. 

Concerning the mean free paths, we have to make difference between the elastic and inelastic scatterings. Under 
usual experimental conditions the elastic mean free path {Id) is smaller than the inelastic one (^in). Here the inelastic 
mean free path is the length of the path an electron travels between two inelastic scatterings {hn=VF'''in)- Mostly 
the important parameter is not /in, but the inelastic diffusive length, that is the average distance an electron can 
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(a) (b) (c) 

FIG. 3: Three point eontact models, (a) An orifice with diameter "d" in an insulating screen between two metallic half-spaces, 
(h) Two bulk regions are connected with a long, narrow conducting neck, (c) Rotational hyperboloid. 

diffuse between two inelastic scatterings: iin=\/^in^ci- If the contact diameter is much smaller than any of the mean 
free paths d <C /ci,^in, we speak about a ballistic contact. In this case the electron travels through the constriction 
without any scattering (except for the reflection on the walls). On the other hand, if d ^ ^ci, the electron makes a 
diffusive motion in the contact, and accordingly we speak about the diffusive regime. At contact diameters exceeding 
the inelastic diffusive length {d » Ljn) the excess energy of the electrons is dissipated inside the constriction, which 
causes a considerable Joule heating in the contact. This limit is called thermal regime. In the following subsections 
these different limits are discussed. In many cases the system is characterized by the Knudsen number, K = d/lei, 
which was first introduced for the problem of the gas outflowing from a tank through a }iole,JSiili^S££i224 but in the 
recent decades it has been used to characterize point contacts as well. 



A. Diffusive regime 

First the diffusive contacts are treated , for which the electric potential (<&) can be determined by classical equations. 
(For a general discussion see lHolml967() If the mean free path of the electrons is much shorter than the dimension of 
the contact (Zci ^ d) then the current density, j is given by Ohm's law in terms of the electric field, E or the electric 
potential, 

j = o-E = -ctV$, (20) 

where a is the conductivity of the metal. Furthermore due to the charge neutrality in metals the continuity equation 
holds: 

Vj - 0. (21) 

If the conductivity is considered to be constant, these two equations result in the Laplace equation for the electric 
potential: 

A$ = 0. (22) 

In this phcnomcnological approach the scattering processes are included in the conductance, which is inversely 
proportional to the mean free path, a ~ This treatment does not distinguish between elastic and inelastic 

electron scattering. In the further treatment of the structures in the dynamical conductance (Sec. lIVp the inelastic 
scatterings play crucial role, therefore, instead of the phenomenological theory the kinetic equations must be used. 
Here the solution of the Laplace equation is presented in order to determine the linear resistance of a diffusive contact, 
and to determine the position dependence of the electric potential. 

Furthermore, we note that the Laplace equation is only applicable if a is constant is space. This assumption is 
certainly not valid if the temperature of the contact neighborhood is inhomogeneous, which happens in the thermal 
regime. This situation is treated in the next subsection. Here we assume that the inelastic diffusive length is much 
larger than the dimension of the contact, thus the temperature of the contact region is constant. 

The Laplace equation should be solved with the boundary conditions: 

^z^±oo)^T^. (23) 
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FIG. 4: The adequate elliptical coordinate system, in which the Laplace equation is solved for an opening like point contact. 



This problem was first solved by MaxwelliMassiii224 Due to the geometry of the problem the Laplace equation is 
most easily solved in an elliptical coordinate system demonstrated in Fig. 01 (see e.g. Morscl953, Simonyil980). The 
transformation between a traditional cylindrical coordinate system (r, z, (p) and the elliptical coordinates (A, /x, if) is 
determined by the following equations of the ellipses and hyperbolas: 

1 and 



a2(l + A2) a2^2 

^2 ^2 



a? (1 — /i^) a^/i^ 



1, (24) 



where a = (i/2 is the radius of the contact; — oo < A < oo; and < /i < 1. Using Eq. (|24|l the old coordinates can be 
expressed in terms of the new ones: 



r = aV(l+^T(l^^Ai^> z = a\^JL. (25) 

In this coordinate system the orifice is given by the A ^ surface, while the isolating layer is defined by /i = 0. The 
next task is to determine the components of the metric tensor in the new coordinates: 

ds^ = dr^ + dz^ + r'^Atp^ 



or, using A and /j,. 



^ a^^l±^A\^ + a^t^^±^dfi^ + a^l + A^) (l - (26) 



The prefactors are the components of the metric tensor in the new (also orthogonal) coordinate system: 

^2 _|_ ^2 ^2 _|_ _)^2 ^ 

gAA = a % , ^2 ' di^t^^a^-. 2"' 5w = + (l ~ M^)- (27) 
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FIG. 5: The potential as a function of the distance measured from the center of the opening type contact along the z axis, in 
the ballistic and the diffusive limit. 



Substituting these into the form of the Laplace operator expressed by the general orthogonal coordinates: 



one can get easily 

A$ = 



1 



V511522533 



a2 (A2 + fi-^) 



d 

dxi 



d_ 



922933 

gii dxi 



d 
dx2 



911933 
922 dx2 



d 
dx3 



d_ 



9$ 



d 



911922 
933 dx3 



9$ 



dfi J d^\{l + A2) (1 - /i2) dip 



Taking the geometry of the problem into account it is obvious that 

- 

dip 



(28) 



(29) 



(30) 



The boundary condition Eq. H23|) can only be satisfied if at A — > ±00 the potential is independent of /x, that is ^ = 0. 
As an Ansatz we extend this condition to arbitrary A values which means that the equipotential surfaces are the 
A = constant surfaces. As these surfaces are orthogonal to the z = plane it is obvious that electric field has no 
normal component at the insulating layer. After this consideration the Laplace equation takes the following simple 
form 



d 

dA 



l + A^ 



d$ 
dA 



0. 



The solution of this equation with the boundary conditions is: 

V 

$(A) = arctan(A). 

TT 

Along the z-axis (/i ~ 1), the potential is changing as 

9[z) = arctan — , 

TT a 

which is plotted in Fig. 

Inside the opening (A = 0) the electric field has only z component: 



^ _ 9$ _ 1 a$ 

^ dz a/i dX 



V 



A=0 «^ VWF 



(31) 



(32) 



(33) 



(34) 
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It should be noted that the field is the largest at the edge of the opening (r — > a) and most of the current flows 
through that region. The total current flowing through the contact is: 



/ = o-y 27rrdr E^{r) = 2aaV, (35) 



and thus the resistance of a diffusive point contact is: 

This calculation was performed for an opening type contact. In case of a channel type geometry with L ^ d the 
potential drops in the channel, and the resistance is easily written using Ohm's law: 

R^^- (37) 
B. Thermal regime 

If the dimension of the contact is larger than the inelastic mean free path, then the electron's excess energy is 
dissipated in the contact region, and Joule heating takes place. In this case the contact center can be considerably 
overheated compared to the bath temperature, and thus the conductivity is also position dependent: a{T{r)). This 
phenomenon can be treated by considering the equations both for the electrical and the heat conduction. The electrical 
and thermal current densities are: 

j = -ctV$, q = kVT, (38) 

where k is the heat conductivity and T is the position dependent temperature. The continuity equations for the 
electrical and thermal current are determined by the charge neutrality and the Joule heating respectively: 

Vj ==0, Vq = -jV$. (39) 

Furthermore, we assume the validity of the the Wiedemann-Franz law: 

^ = Ca, (40) 

where £ — {■nks)'^ /ie^ is the Lorentz number. Combining these equations we get: 

V(ctV$) = (41) 
^V(CTVr2) +cr(V$)2 = 0. (42) 
From these equations it can be easily seen that the temperature is generally related to the electric potential as: 

r2(r)= const. (43) 

The constant term is determined by the boundary condition: 

T(0^±oo)^rbath, (44) 

thus: 

T\v)^TU + ^-'^- (45) 
In an opening type contact the potential is zero at the contact surface, so the local temperature is written as: 

Tpc = ?llth + 2, (46) 

According to this relation a bias voltage of 100 mV already heats up a contact from cryogenic temperatures (~ 4K) to 
room temperature. This feature has crucial importance in experimental studies. If the measurements are performed on 
a large contact or at high bias voltages, it can easily happen that the I{V) curve presents the temperature dependence 
of the conductance instead of the spectroscopic featuresr^SS^SiiiiSS which will be discussed in Sec. IIVI 
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C. Ballistic regime 



If the linear dimension of the contact is much smaller than the minimum of the mean free paths the electrons are 
travelling along straight trajectories, and in the first approximation no scatterings are taken into consideration. This 
system is described by the equations of semiclassical dynamics, so the space and momentum dependent distribution 
function, /p(r) is to be determined beside the electric potential, $(r). In this limit the electrons "remember" which 
side of the contact they are coming from, thus the distribution function at position r can be expressed as a sum of the 
terms corresponding to electrons coming from the left or the right side of the contact. The distribution function 
can be determined by solving the homogeneous Boltzmann equation which contains no collision integral(^^iiiiJii22i 



^ dr dp 



eB^^q^ = 0, (47) 



where E is the electronic field, E = —^^7^, and Vp is the velocity of the electron with momentum p. Note that the 
electronic charge e is positive, so the charge of an electron is — e. Far from the constriction the boundary condition 
for the distribution function is that it has to tend to the equilibrium distribution function. 

hm /p(r) = /o(£p), (48) 

|r| — >oo 

where fo{£p) — (e^r^ + 1) 

If we apply V voltage on the sample, the boundary condition for the potential $ takes the form 

$(z^±«)) = Ty. (49) 

The solution of the Boltzmann equation can be found by using the trajectory method. "^iSSiJaiiSiiSiiiSZi (The trajectories 
of the electrons are considered to be straight because we are interested in the limit e|y j/ep ^ 1 and the current we are 
searching for, depends linearly on the electrical field. It can be shown that the bending of trajectory only contributes 
to the current in the higher order of the field.) In this way 

/p(r) = /o ^£p + e j Edlj (50) 

holds combined with the charge neutrality condition 

2eyd3p[/p(r)-/o(£p)] =0, (51) 

r 

where J denotes the integration over the trajectory of the electron coming from the distant reservoir to the contact 

—00 

region at point r. The distribution fmiction of an electron satisfying the conditions H51|l . (|5U|I and 148|) takes the form 
in the linear order of the field 

/p(r) = /o (^ep - e$(r) + ^r^p, r)) , (52) 

which is constant along the trajectory as the energy Sp — e<I>(r) is conserved. The expression eVri{p,r)/2 describes 
the different value of the chemical potential in the right and left hand sides of the constriction {fi = mIv=o ± eV/2): 

I s _ J +1 for electrons arriving at r from the left reservoir 
^ 1—1 for electrons arriving at r from the right reservoir 

(For an illustration see Fig. 0) That can be also expressed in terms of the solid angle, ri(r) at which the opening of 
the contact is seen from point r: 

'71P, r j - I _sign(z) , if -p ^ n{v) ' 
A simple visualization of the distribution function is presented in Fig. 
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FIG. 6: Different electron trajectories arriving at point r. The values ofrj{p,r) — ±1 are also indicated. 



0=V/2 




= -V/2 



z 

-► 



FIG. 7: Visualization of the momentum distribution function at three different points in the case of ballistic limit. The shadowed 
areas represent the occupied states at function of the direction of the momentum, p. The dashed circle represents the equilibrium 
Fermi surface in case of completely isotropic distribution. The volume of the shadowed areas are independent of the position r 
(charge neutrality). 



The potential can be derived by substituting the form H52() of the distribution function into the neutrahty condition 



2e / d-^p 



fo[£p- e$(r) + ^viP, r) ) - /o(£p) 



0, 



(55) 



where the integral for equienergetic surfaces is easily expressed by the solid angle ri(r), and just an energy integral 
remains: 



de 



0, 



(56) 



where the upper and lower signs correspond to the cases z > and z < 0, respectively. The equation can be solved 
easily by using the identity 

(57) 



(58) 



(/o(e + a) - /o(e))de = -a 
for an arbitrary energy shift a. We can get the potential in the whole space as£iagijaii£ii2Jiaffi 



$(r) 



V 
"2 



2n 



sign(2;). 
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Here fi(r) is the solid angle at which the contact is seen from the position r. The potential depending on the distance 
measured from the contact along the z axis is shown in Fig. [S] (solid curve). 
The current density, j is determined by the distribution function as follows: 

j(r) = -2e 1 ^vp/p(r). (59) 

The current through the contact is calculated by integrating the z-component of the current density over the area of 
the contact, A: 



1 = J dAj,(z = 0), (60) 



At the contact surface half of the electrons go to the left with a distribution function /o(e — eV/2) and half of the 
electrons goes to the right with a distribution function /o(e + eV/2), thus the current is written as 



where the integral J dS^ is taken over an equienergetic surface in the p-space. At low enough temperature and voltage 



{k^T <C Ef and eV <C £f) the expression J (/o(£ + ^) — /o(e ~ ^))de equals eV according to the identity (|S7|), so 



the current is written as: 

~ ?i(27r)3 



cos(79))fs • V, (62) 



where A is the area of the contact, ^f is the area of the Fermi surface, cosi? = Vz/vp and (...)fs — J d5'F(...)/<S'F 
denotes the average taken over the Fermi surface. The conductance of a ballistic point-contact given by the formula 
(I62|l is known as the Sharvin conductanceiSiaaaiiiSSi For free electron gas the Sharvin conductance is simplified as: 



D. Intermediate case between the ballistic and diffusive limit 



In the intermediate region between the diffusive and ballistic regime an interpolating formula can be set up by 
numerically solving the Boltzmann equation for arbitrary ratio of the contact diameter and mean free path, ;|Wexieri966 

R = l/d-^+Til/d)^, (64) 
oiraa ad 

where ^{l/d) is a numerically determinable monotonous function, r(0) = 1; r(oo) = 0.694. Note that the first term 
is exactly the Sharvin resistance by putting the Drude resistivity into the formula, p — mvp /le^n, thus it is actually 
independent of /. This formula provides the possibility to estimate the contact diameter from the contact resistance. 



IV. SCATTERING ON SLOW TLSS IN POINT CONTACTS. 



The voltage applied on a point contact results in a nonequilibrium distribution of the conduction electrons in 
the contact region. An electron coming from the right reservoir has an energy larger by eV than those coming 
from the left reservoir. This energy can be released through inelastic scattering processes, which can happen in 
such a way that an electron that has already crossed the contact is scattered back through the opening. This so- 
called backscattering reduces the current, thus the energy dependence of the scattering probability can be traced 
by measuring the nonlinearity of t he current as the bias voltage is varied. This phenomenon was first used by Igor 
Yanson^iS§sai22i (for a review see lYansonl986al) to study the phonon spectra and the electron-phonon interaction 
and since that it is widely applied. 

In this section we discuss the nonlinearities in the current voltage characteristics due to the scattering on slow 
TLSs, which show strong similarities to localized phonons. There is, however, an essential difference between the 
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two cases, while the phonons contribute to the voltage region up to 30 meV, the TLSs manifest themselves below or 
even well below 1 meV due to their characteristic energies. The main contribution to the I{V) characteristics comes 
from the close region of the contact, therefore a microscopic process like a transition of one atom from one position 
to another in the contact region occurs as a significant measurable change in the current. Considering the TLSs the 
main advantage of the point contact spectroscopy is that in case of small enough contacts one can investigate just 
one or few scattering centers. 

The current through the contact can be derived by solving the classical stationary Boltzmann equation 

where /(p, r) is the collision integral for the electron with momentum p at position r, which is assumed to be local 
in real space. 

In the case of bulk phonons the collision integral is 

/ph(p, ^) = E / I [/p + q(l - /p)(^q + 1) " fpi^ " /p+q)^q] '^C^p + q " £p " + 

[/p_q(l - U)N^ - /p(l - /p-q)(iV^ + 1)] <5(ep-q - - UJ^)^ , (66) 

where Wq is proportional to the squared matrix element of the electron-phonon coupling, is the phonon occupation 
number for momentum q and energy Wq in phonon branch a. There is another term of the collision integral which 
is due to the elastic (impurity) scattering but that is not given here. The collision integral vanishes far from the 
point contact because well inside the electrodes thermal equilibrium exists. The situation is further simplified when 
the following assumptions are valid: (i) the electron-phonon interaction is homogeneous in the real space; (ii) the 
phonon distribution corresponds to thermal equilibrium. In that case the information on the contact geometry and 
also the elastic scattering due to impurities in the dirty limit can be expressed by a geometrical factor K{p,p') in 
the expression of the current which depends only on th e momenta p and of the incoming and outgoing electrons. 
That factor is frequently called as the i^T-factor (see e.g. lKulikl98llAshrafl98^ . In the case of phonons the K-ia.ctor 
for a ballistic opening like contact is written as: 



\p'zP- 



KiP, P') - Tir^^^J^\^(-PA), (67) 



where p p' are the unit vector parallel with p and p', while Pz and p'^ are the z components of these unit vectors, 
where the z direction is the axis of the contact. The A'- factor can be also calculated for a diffusive contact: 

K{p,p')^^^^-f{pl^pf), (68) 

where Id is the elastic mean free path, and d is the contact diameter. 

At zero temperature the logarithmic derivative of the resistance can be expressed by a function G{uj) = a'^ (uj) F (uj) , 
where F{uj) is the density of states of phonons and a{uj) is the modified electron-phonon coupling which contains the 
if-factor. Without that factor it is just the electron-phonon coupling strength and in that case G{uj) is the Eliashberg 
function known in the theory of superconductivity. The main result is 

In case of single crystals G{lu) depends also on the orientation of the crystal, which makes it possible to study the 
anisotropy of the phonon spectrum. Measuring the phonon spectra were determined for many different 
cases and the structures in them were identified as the details of the spectraJ^ii2ii2Sl£iii22^ Those structures, however, 
are superimposed on a continuous background. 

Until now, it is assumed that the phonons are in thermal equilibrium, which is the correct assumption as far 
as the phonons at the contact region arrive from a large distance without any collision where the thermal bulk 
distribution is realized. That is, however, not the case when in the contact region the phonon mean free path is 
comparable with the contact size. Then the phonons generated by the non-equilibrium electrons in the contact 
region cannot escape thus the phonons are also out of equilibrium. The measure of non-equilibrium is limited by 
the energy relaxation time of the phonons Tph. The non-equilibrium localized phonons contribute to the background 
resistancei^™^'^''^'^''^^'^^i^'^'^Q''^'^'^''^^^°iJ™"™^^^° That can be studied by measuring the amplitude of the background as a 
function of the frequency of the applied vohago^""'^^^^'^!^'^'^'^"'^^^^'^!^'^"^'^"^''^^^^^ and at high frequency (wTph > 1) that 
is decreasing. These phenomena will be studied in the context of scattering on TLSs as well. 
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A. Scattering on TLSs in a ballistic contact 



In the following the Boltzmann equation will be solved considering only the electron-TLS interaction. It will be 
assumed that the density of the TLSs is low thus an electron scatters only once in the orifice region and the contribution 
of the different TLSs are additive. First the ballistic limit is treated, where both the elastic and inelastic mean free 
paths are larger than the size of the contact. 

The collision term related to a single TLS is very similar to the case of phonons, the momentum is, however, not 
conserved. Furthermore, the TLS has only two states with population n+ and n_ + n_ = 1) corresponding to 
the TLS energy levels i?+ and E- (see Sec. The situation is different from those cases of phonons, where the 
phonon mean free path is large compared to the size of orifice, Zph 3> d and the phonon distribution is the bulk one 
which is in thermal equilibrium with the contact regions. In the case of TLSs the interaction is localized in space, 
thus the TLSs in the contact region are decoupled from the bulk region and they can be considerably out of thermal 
equilibrium due to the voltage drop at the contact. Therefore, the situation is similar to the special case of localized 
phonons with /ph < d shortly discussed in the introductory remarks of this section. 

The collision term similar to Eq. (|6til) can be written as the sum of inelastic and elastic terms as /R,TLs(P;r) = 
/R_TLS,in(Pj r) + /r,tls,oi(Pj t). In the inelastic case the TLS jumps to its other state due to the scattering, and thus 
the energy of the electron changes hy E = E^ — E^ . The strength of the interaction is characterized by Wpp' , which 
is obtained from the interaction matrix elements (Eq. I12() by Fermi's golden rule as 

Wpp' = Y \{P'E+\H,-TLs\p\E^)f = ^ |(p,i?_|i/e-TLs|p/i?+)|' - ^ \2f^^V^p,\^ . (70) 

The interaction strength is symmetric in the incoming and outgoing momentum, and in the neighborhood of the 
Fermi surface it depends only on the unit vectors p and p', i.e. W^pp' — Wpp' — Wp'p. The collision integral for the 
inelastic scattering is written as: 



/R.TLS,i„(p,r) ^ 6{r-K) J -0^W, 



|/p'(l - .fp)n-Siep- -ep-E)+ /p.(l - /p)n+5(ep- - Sp + E) 

-/p(l - /p')"-'^(£p - V /p(l - fp')n+S{Sp - V + i?)|, (71) 

where R is the place of the TLS. 

The current correction due to the elastic scattering processes can be written similarly with two remarks: (i) in an 
elastic process the energy of the TLS does not change, thus E — must be inserted; (ii) the scattering cross section 
is different for the TLS being in the "+" and the "— " state. According to the matrix elements in Eq. H13|l. the two 
scattering strengths are given as: 

W±, = ^|(p,i?+|i/,_TLs|p',i?+>l'-^I^PP' + (^'~"')^PP'l' 

W--, = ^|(p,i5_|i/^_TLs|p',i?-)l' = YlC'-(/^'-'^')^PP'r (72) 
After introducing these notations the correction to the collision integral for the elastic scattering is: 

/r,tls,ci(p, r) c Sir -K)J {w±,n+ + Wpp,n^) [/p.(l - /p) - /p(l - /pO] 5(ep - v)- (73) 

The electron distribution function /p(r) and the electric potential is expanded in terms of d/im- At a large distance 
measured from the orifice the potential is constant and lim|r|^oo '^(r) — ±V/2. Thus, at |r| ^ li^ the electrons are in 
thermal equilibrium. The distribution function /p(r) and the potential <I>(r) can be expanded as 

/p(r) = /(")(r) + /W(r) + ... 

$(r) = $("^(r) + $(i)(r) + ..., (74) 

and similarly the electric field 

E(r) =E(o)(i-)+E(i)(r) + ..., (75) 
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where the upper indices (0), (1) label the order in the strength of the electron-TLS coupling. The zero order terms 
have been previously calculated and the following treatment is restricted to the first order terms. The previous results 
are 

(r) = /o f ep - e$(r) + ^rj{p, r)] , (76) 



and according to Eq. 



$(0)(r) = -^ fl-^)sign(z), (77) 



2 V 27r 



for a round orifice perpendicular to the z-axis, where /o is the Fermi function. The first order terms in the Boltzmann 
equation given by Eq. (|65|l can be written as 

- eE™,„^^ = eE« H- .^(P.,). (78, 

where the label (0) of the collision term indicates that the collisions are calculated with the distribution functions of 
zeroth order, /(°^. 

The change in the potential is determined again by the neutrality condition given by Eq. H51|l for f^^\ which is 



'^'P /W(r) = 0. (79) 



(27rfi) 



The boundary conditions are 



and 



hm /(i)(r) = 0, (80) 

\r\~^OQ 



lim $(i)(r) = 0. (81) 

|r| — ^CJO 

In order to determine /*^^^ and the trajectory method is used (see e.g. lJansenl980() dealing with phonons. Electron 
trajectories are considered where the electron comes from the left or right and arrives to the plane of the orifice at 
time r = at point r with momentum p to calculate the current flowing through the orifice. These trajectories 
correspond to zeroth order. Moving along the trajectory r(r), p(t) the distribution function varies as 

9/(i)(p(r),r(r)) _ 0p(r)9/W(p(r),r(r)) ^ 9r(r) a/(i)(p(r), r(r)) ^ 
dr dr 9p(t) dr dr^r) ' 

where ^^p- = — eE, and = Vp. This equation can be integrated and using the Boltzmann equation one finds 

for /(I) at r = 



Here the collision term describes the adding or taking off electrons to the trajectory arriving at the contact at r = 0. 
In the expression of f^^^^ only the terms linear in the voltage are kept, thus — — g'p^^'^^^^^ in the above equation can 
be approximated by the zeroth order term Vp^j^.. Taking into account that dr = Vpdr the first term of the integral 
can be performed, thus 



/«(p,r) = -e$W(r)|^+ J dr /(°,\(p(r), r(r)). (84) 
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The equation of electrical neutrality taken at r = with p(t = 0) = p and r(r = 0) = r combined with the equation 
above determines $(^)(r) 

I A I dr4"i(p(r),r(r)) 



^*^^Hr) = "7 d3p of^ ■ (85) 



The general expression of the electrical current flowing through the surface of the contact is 

^ = -2e/d^^^/^«^/(p,r), (86) 

where d^g is the surface element of the contact and the integral is taken over the orifice, the factor 2 is due to the 
electron spin. The change SI in the total current due to the presence of the TLS can be separated as 

I{V)^Io{V) + 5I{V), 

and SI{V) can be further split, whether the electron scattering is elastic or inelastic 

Making use of the expression H84|) for /^^^ only the collision term contributes to the current because the other term 
is even in the momentum. The integral according to the time can be transformed to the one along the path using 
dr = ^ where ds is the element of the path. The expression for the change of the current due to the presence of 
collision on TLS is 

— oo 

The path of the integral is changed due to the collisions but it may contain reflections by the insulating surface of the 
contact. The next step is to introduce the collision term due to a single TLS at position R, given by Eq. (|71|l . As the 
paths arriving at the surface element of the orifice d^ g are straight lines, the volume element for the scattering event 
is d'^r = d^gdscos^ where 6 is the angle between the path and z direction. Now the integration over different paths 
and position of the collision gives 

SlR = -2e J ^i!^^cos0./R,TLs(p) (88) 

where /r,tls is defined by jR_TLs(Pjr) = 5{r — R)/r,tls(p) and the momentum integral is restricted to the sohd 
angle f^R at which the contact can be seen from the TLS at position R. Furthermore, fp/fp = —cos9 sign(i?^) as 
the electron passes the contact from the TLS (for an illustration see Fig. IS)). The final expression is obtained using 

f^2^nj3 = f'o^^^de where dflp is the solid angle element in the momentum space and is the conduction electron 
density of states for one spin direction. Assuming a spherical Fermi surface 



(5/r 2e sign(i?^)eo J ^ de cos^ 6* • /r,tls (p 



)• (89) 



peOp 



1. Current correction related to inelastic scattering 



Using the expression of the coUisison term /R,TLS.in (I71|l the change in the current due to the inelastic scattering 
can be written as 

^Ran = -^sign(i?2) / j de cos^e / -^;^Wpp' / de 

PGOr 

{ [/°,(1 - /^)n_ - /^(l - f^,)n+] 6(8 + e') + [/^,(1 - /^)n+ - /°(1 - f^,)n.] 6{e - E - e')}, (90) 
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TLS 



FIG. 8: The orifice with the path of the electron scattered by the TLS on the right. The current flows through the surface 
element d? g and the angle 6 is also indicated. 



where the distribution functions f'^ are given by Eq. (|76|l and taken at r = R, and the dimcnsionless notation 

Wpp' = hglWpp' 



(91) 



is introduced. The potential *&(R) appearing in these expressions can be dropped as the energy integral variables can 
be shifted. The essential contribution is given by ^ri{p, R) = which has a positive sign if the electron at r = R 

arrives from the left contact and a negative sign if it is arriving from the right contact. 
The final expression is 



Si- 



ll,) 



2e 
IT 



sign(i?^) 



47r 



-Wt,i 



ds cos2 e I ^ I de' 



foie' + ^ry(p'))(l - /o(£ + ^vip)))n- - fo{e + ^^(p))(l - foie' + ^vip')))n+ 5{e + E - e') 



foie' + ^r;(p'))(l - foie + ^r;(p)))n+ - Us + ^vipW ~ fo{e' + ^^(p')))"- S{e - E - e') 



(92) 



This expression can be divided into two parts. If the electron with momentum p' arrives from the same side of the 
contact (|r| oo) as the unscattered electron with momentum p, thus '7(p)f7(p') = 1, then the electron is scattered 
forward, while for ?7(p)»7(p') = —1 it is scattered backward (see Fig. I^)). The total current can be divided according 
to that, thus 



5Ib 



R,in — ^/r, in, for + 0-'R,in,back 



(93) 



The forward scattering cancels out. That can be seen by looking e.g. the first and the last term in the previous 
expression and changing e ^ e'. 



5h 



R,, in, for 



= 0. 



(94) 



In the calculation of the backscattering the case sign(i?2) = 1 is taken first. Then the electron contributing to the 
current originally comes from the left (»/(p') = 1) and after the scattering goes from right to the left (?7(p) = —1). 



«R.i 



R,. in. back 



2e 



£ / de cos^ 



d^p, 



de' Wpp> 



-P 'efiR 
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|{/o(e' + ^)(1 - Me - - fo{s - ^)(1 - fo{s' + ^))n+}% + - e') 

+ {Me' + ^)(1 - /o(£ - - /o(£ - - /o(e' + ^))n-}s{e -E- (95) 

The factor due to the angular integrals J cos^ 9 J -^^wppi plays the role of the geometrical factor 

for the bulk phonons. Only the backscattering contributes like in the case of phonons as it can be seen from Eq. H67() 
for ballistic regions where the if-factor contains the factor Q{—pzp'^). 

In the following the geometrical factor is introduced which depends on the position of the TLS, R, and it is 
independent of the strength of the interaction 



i / ^cos^. / ^ 
w J At: J An 



^1^ = - / -rrcos^^ / ^^pp', (96) 



where 



If the momentum dependence of the interaction is ignored then Wppi = w and 

... / '-^r,o.H I PS, 

pgOr -p'eJ2R 

is a pure geometrical factor. Further on this simplification is assumed. If the TLS is positioned on the geometrical 
axis of a circular contact then the above integral can be evaluated as 

In the second part of the integral in Eq. (|95|) the variables are changed as e ^ e' and the integral with respect e' 
is performed, furthermore, in some of the integrals the variable is shifted The result is 

«Ran,back = ^^^R J /o(e) { [/o (^ " ^ + eV) - Me - E - eV)] n_ 

+ [Me + E + eV) ~ Me + E - eV)]n+y (100) 
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Now the TLS on the left is considered (sign(i?2) = —1). Then for backscattering the signs of ?7(p) and ri(p') are 
reversed (see Fig. |5Jd) That is equivalent of changing the sign of V {V ^ — F) in Eq. (|95|l and then the expression 
in the curly bracket also changes its sign as it can be seen by inserting e ^ e'. That negative sign is cancelled by 
sign(i?2) = — 1, thus the total current is unchanged. Therefore, the current does not depend on whether the TLS is 
on the right or left as it should be. 

After inserting the Fermi functions and performing the energy integral, the following result is obtained for the 
inelastic current correction: 



(101) 



where 



2e \eV + E 



coth( 



eV + E 
2kBT ^ 



2e 

D = -j-wKj^ 



eV -E 



coth(— — — ) - 1 



2kBT 



These expressions are simplified at T = as 



E -eV 



-E 



2 

eV 



coth( 
coth( 



E-eV ^ 
2kBT 

-E - eV 



2kBT 



e\V\ < E 



2e I '2e\V\n+ 

Sln.in = -^wKu.sign{V) 

' 2En+ + {e\V\ ~ E) e\V\ > E 



(102) 



(103) 



If the TLS was in thermal equilibrium with the bath, at zero temperature n+ = and n_ = 1 would hold. In that 
case the voltage dependence of the conductance is easily evaluated, it contains a jump-like decrease at the excitation 
energy of the TLS: 



dG 



R,in 



R,i 



dV 



^^'^wK^-Q{e\V\~E), 



(104) 



where Q{x) = | J 5^^". That result means if the TLS is in the ground state n_ = 1 then the TLS can be excited 

only by a voltage e\V\ > E, where the backscattering reduces the current. The second derivative has a Dirac delta 
peak at the excitation energy 



d^Sh 



R,i 



dSG 



R,in 



2e^ 



dV 



= —w 



Krisign{V)-6{e\V\-E). 



(105) 



This result is similar to the phonon result, that is the shows the structure of the excitation spectrum of a scatterer 
in the contact region. However, for a TLS situated in the contact region the occupation of the upper level is not zero, 
as the system is driven out of equilibrium. Due to this feature a nontrivial contribution occurs in the background 
resistance beside the spectroscopic feature at the excitation energy. This can be evaluated after calculating the voltage 
dependence of the occupation number, n^{V), which is performed later in this section. 



2. Elastic scattering 

A similar calculation can be performed for the current correction due to the elastic scattering on TLSs by using 
the collision integral /r,tls,oi(p, r) (|75|l instead of /R,TLS,in(p, r) EJ. The formulas for the elastic scattering can 
be derived easily by inserting E — and modifying the scattering strength in the results for the inelastic case. For 
elastic scattering the strength of the interaction can be different for the TLS being in the "+" and "— " state, thus in 
Eq. (|102|l wi^R must be replaced by [w/Cr.]" and [wXr]^ in the expression for C and D, respectively. After these 
considerations the result is: 

SIr^cI = (n- • [wKr]- + n+ ■ [wKr]+) . (106) 

Note, that in the elastic case the voltage dependent "coth(- • •)" terms cancel out, thus the only voltage dependence 
comes from the voltage dependence of the occupation number. If the TLS is in thermal equilibrium with the bath, 
then n+ = 1 — n_ = holds, thus the elastic current correction is a voltage independent constant term. Similarly, if 
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the scattering cross sections for the two states of the TLS are equal, then the voltage dependence cancels out due to 
the condition + 7i_ = 1. In the following we use a simplified notation for the elastic current correction: 

(5/ei = -l/(7"n_+7+n+), (107) 

where 7+ and 7" are considered as the reduction of the conductance due to the TLS being in the "+" and the "— " 
state, respectively. K™">^1986 rpj^g 

average change in the conductance due to the elastic scattering is characterized by 
t"'^ , which is written as follows according to formulas (|106l 1981 1911 1721) : 



7+ + 7" _ 2e^ 
2 ~ ~ 



(27r)2i^Rf^^([y°]%[(M^-z.2)F^]^), (108) 



where and are the the interaction matrix elements from Eq. (|ll|l assuming isotropic scattering. The difference 
between the conductances corresponding to the two states of the TLS are written as: 

2e2 

7+ - 7" = —(2T,fK^Ql(\^^ - 4i.2)^o^._ (^09) 

Again, the voltage dependence of the current correction can only be determined after calculating the occupation 
number as a function of bias voltage, which is done in the following part. It is important to note that in case T^" = 
the scattering strength is symmetric for the two states, thus 7"*" = 7". 



3. Calculation of the occupation number 



In the following n+ will be calculated in the nonequilibrium case, where it is determined by the interaction with 
the nonequilibrium electrons. The transition probability from the "— " to the "+" state due to an electron scattered 
from momentum state p to p' has been considered in the collision term H71|l and it is: 



W+-ip'p) = T4^p.pn_/p(R)(l - /p'(R))<5(£p - e^, - E), 

and similarly: 

W^R + (P'P) = W^p'p"+/p(R)(l - /p'(R))'5(£p - ep,+E). 
The transition probabilities for the TLS are obtained after integrating for the electron momenta: 

W^iC" ="T^/^ / ^^P'P / de{/o(e-e$(R) + ^,Kp))(l-/o(eTi?-e<f(R) + ^ry(p'))}, 



(110) 



(111) 



where the expression ((7^ for /p(R) is used. 
The kinetic equation for n+ is: 



dn+ 



R 



W^R + - 



(112) 



(113) 



In the Fermi functions the electron momentum is only involved in ?7(p), which tells which reservoir the electron at 
position R with momentum p is coming from. Therefore the integral can be performed by separating the regions of 
integration into four cases with respect to the possible values of T]{p) and ''7(p')- In all of these four regions the factor 
containing the Fermi functions is independent of Qp and Q'p. There is an essential simplification if wppi is replaced 
by a momentum independent averaged one, w. Take e.g. the case where the TLS is on the right hand side of the 
contact {R^ > 0) and introduce the solid angle 17r at which the opening can be seen from point R. In this case the 
integrals for the different regions are written as: 



(•••)d£ 







if 


'7(P) 


= 1, r;(p')-l 


fl - ^ 




if 


'7(P) 


= -1, v{p') = -l 


^fl - 


47r > 


if 


v{p) 


= 1, r;(p') = -l 


ilRfl _ 

4tt V-^ 


HrN 
47r 


if 


vip) 


= -1, ?7(P')-1 



(114) 
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This simplification is justified only after taking an average over large number of TLSs but that simplifies the calculation 
and makes the result more transparent. In this simplified case the kinetic equation for the TLS is obtained after 
calculating the energy integrals for the different cases. (As before, e$(R) is eliminated by shifting the integral 
variables.) The final result is independent whether the TLS is on the right or left, and it can be written in the form: 



^ n-A - n+B = A-n+{A + B), 



(115) 



where the coefficients A and B are: 



eV + E { ^ ,eV + E, 
coth(- 



eV + E 



coth(- 



^ + \^ coth( 




2kBT ' 



E 



B = 



cQth( ) - 1 



2kBT 



coth( 



-E 



2kBT' 



(116) 



respectively. The relaxations are of Korringa types, which depend on V only if the two electrons involved have different 
chareacters in terms of ?7(p). The stationary value of n+ is obtained as 



A 



^ A + B' 

which is an even function of V^. At T = this expression simplifies essentially: 



(117) 



e\V\ 
e\V\ 



< E 
> E 



(118) 



where k = ^(1 — is a continuous function of V. Far from the opening k ^ thus n _)_--> at voltages for 

which {e\V\ — E)/E <C as the electron gas is in thermal equilibrium. For k = const, and {e\V\ — E)/E 3> k,^^ 
^ i as for large enough voltage the levels are equally occupied = n+ = j. If the TLS is in the middle of the 
contact, then k — 1/4. 

If the TLS is far away from the contact then it is in thermal equilibrium, which can be obtained by taking JIr — > 
and the deviation from equilibrium is large when the TLS is in the middle of the contact region JIr — 27r. Here other 
relaxation mechanism different from the scattering of electrons is not taken into account as the generation of bulk 
phonons is very weak as the relevant phase space is very small at low energies. At large concentration of TLSs the 
TLSs are interacting and the collective effects may modify the stationary values of the occupation numbers. 



4-. Conductance with inelastic and elastic scattering 

Now the inelastic and clastic contributions to the current and the conductance are calculated using the stationary 
occupation numbers obtained. The expressions will be given when the TLS is just in the middle of the round opening, 
thus r^R, ~ 2tt {k — i), but the calculation for arbitrary can be easily performed as well. 

At zero temperature the inelastic current correction is obtained by inserting the stationary value of given by 
Eq. ifTTH)) into (jTU^ : 

2e . ^ f e\V\ <E ^ ^ 

«R=o.i„ = --wK^=osigniV) I ^1^1 _ ^1^1 ^ ^ , (119) 

which is a continuous function of V. After differentiation with respect to the voltage the correction to the conductance 
is obtained as: 

9<^/R^o.in Jo e\V\<E 

= = -ir"''^-° I 1 + 2(^^ e\V\>E^ (120) 
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and the second derivative of the current is 



dV 



e\V\)-Qie\V\-E) 



AE' 



{E + e\V\f 



(121) 



At positive voltages the second derivative of the current shows a negative Dirac delta peak at eV = E, which 
reflects the energy spectrum of a single TLS. Above the excitation energy a positive background is seen due to the 
nonequilibrium distribution of the TLS occupation number. 

As a next step the elastic contribution is determined at T = for a single TLS positioned in the middle of the 
contact. The change in the conductance due to the elastic scattering can be calculated by differentiating Eq. (|107|l . 
The elastic current correction contains a linear term, which causes a constant, voltage independent reduction of the 
conductance. Experimentally it is hard to separate this constant contribution, thus we calculate only the voltage 
dependent part by subtracting the zero bias conductance: 



SG 



dSL 



R=0,cl 



dSL 



R=0,cl 



R=0,cl 



dV 



dV 



-(7+ -7") 





1 

2 



e\V\ 
e\V\ 



< E 
> E 



The second derivative of the I{V) curve is 



dV 



1 

= -e(7+ - 7-)sign(y) { -6{E - e\V\) + Q{e\V\ - E)^^^-^ 



(122) 



(123) 



Again, a Dirac delta peak reflects the spectrum of the TLS, and a continuous background arises at e\V\ > E due to 
the nonequilibrium distribution. Contrary to the inelastic case, the Dirac delta peak and the background have the 
same sign; furthermore, in the elastic case the sign of the peak can either be positive or negative depend ing on th e 
sign of 7"*" — 7~. The amplitude 7+ — 7" is related to the universal conductance fluctuation 2e^//i (see e.g. lLeel986l) . 

For arbitrary position of the TLS the occupation number is zero at = and 1/2 at ^ 00. Therefore, the 
total amplitudes for the change in the conductance in the elastic and inelastic case can be generally written as: 



AG 



1 



R.cl 



(7+ -7-) 



2e2 



(124) 
(125) 



where isotropic scattering was assumed, and the formulas H7()ll91ll?^ll()9|l were used. The results are given in the unit 
of the universal conductance quantum, 2e^/h. According to Eq. Q for a highly asymmetric TLS, where the energy 
splitting A is much larger than the transition term Aq, the equations /i ~ 1 and ~ hold. In this case the inelastic 
term is suppressed. In the opposite case, where A <C Ag the elastic term is suppressed as /i ~ ~ 1/2. It must be 
also noted that the inelastic term depends only on the matrix element , whereas the elastic term is influenced both 
by and . The inelastic term can be roughly estimated as AGR^n ^ 0.02 2|- using K^^o = 1/6, = v = 1/2, 
qqV^ 0.1 (see Vladarl983a Vladarl983c). The amplitude sharply drops by moving the TLS further from the orifice 
than its diameter because of the geometrical factor. The estimation of the elastic term is more difficult as is 
strongly model dependent. 

In the following the effect of a TLS is considered at finite temperature. The formulas for the conductance and the 
second derivative of the I{V) curve can be explicitly calculated using equations IjlUll I1031 11061 [TTCl I117|) : however, 
these equations are very complicated, thus the results are demonstrated by figures. Figure 1101 shows the voltage 
dependence of the conductance for the inelastic and elastic case respectively. In Fig. the second derivative of the 
I(V) curve, the so-called point contact spectrum is presented at various temperatures both for the elastic and inelastic 
case. In Fig. E| both the inelastic and elastic contributions are compared for a TLS positioned in the contact center 
(r^R — 2-k) and for a TLS being farther away (f^R — 2-k jh). 



B. Slow TLS in a diffusive contact 



As it has been discussed in Sec. lIII Al the phenomenological treatment must be replaced by a theory based on kinetic 
equation, where effect of the inelastic scatterings are taken into account in the distribution function. In the following 
the elastic and inelastic scatterings are treated in equal footing. As a first step the distribution function is determined 
in the presence of elastic scattering and next the contributions of the TLS are treated as a perturbation. 
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eV/E eV/E 

FIG. 10: The contribution of the elastic and inelastic scattering to the resistivity in the ballistic limit when the TLS is in the 
center of the contact (SGn.^o,in and 5Gr=o,c1, respectively). The temperature is kT = O.l-B. The sign of the elastic contribution 
can be either positive or negative depending on the sign of (7^ —^~). On the other hand, the inelastic scattering always causes 
a decrease in the conductance. 




FIG. 11: The second derivative of the I{V) curve for inelastic and elastic scattering at various temperatures for a ballistic point 
contact (dSGn.^o.in/dV , and 95Gr=o,ci/9V, respectively). The inelastic contribution shows a negative peak at the excitation 
energy of the TLS and a positive background tail at higher voltages. The elastic contributions similarly shows a peak at eV = E 
and a background tail at higher voltages, but here the sign of the peak and the tail is the same. In the elastic case the sign of 
the signal can either be positive or negative depending on the sign of (7^ — ). 



Consider first the elastic scattering in the limit Id 0, where the resistance is very large. Because of the very strong 
elastic scattering, the electrons are immediately redistributed concerning the direction of their momenta. Therefore, 
in the limit lc\ ^ the distribution function depends only on the energy of the electrons, /p(r) — > /e(r). As the 
electron arriving at the contact is either coming from the left or the right reservoirs and the electron energies are 
changed due to the external potential <&(r) the distribution function must be the superposition of the distribution of 
the electrons coming from the left or right with amplitudes ao(r) and 1 — Q;o(r), respectively: 

/,(r) = ao(r)/o(£p - e$(r) + f^) + (l _ ao(r)) /o(ep - e^r) - ^). (126) 

Using the charge neutrality condition (J ff{r)de — J fo{ep)de) ao{r) can be determined by a similar treatment used 
earlier in case of Eq. (|58|) : 



(127) 
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FIG. 12: The second derivative of the liV) curve for inelastic and elastic scattering for two different positions of the TLS: in 
one case the TLS is in the contact center, in the other case it is farther away on the contact axis, so that the opening is seen 
in a solid angle f^R = An/IQ. The TLS being farther away from the contact is closer to thermal equilibrium, which is reflected 
by the reduced amplitude of the background tail with respect to the spectroscopic peak in the inelastic contribution. 



CO 




FIG. 13: Visualization of the distribution function in the diffusive limit. Without inelastic processes the smearing of the two 
steps are due to the finite temperature, T . 



This formula is valid in case of arbitrary geometry because only the shape of the potential function contains information 
about the details of the geometry. For an opening type circular contact the potential was determined in Sec. lIII-^ bv 
solving the Laplace equation using a hyperbolic coordinate system. 

The distribution function contains two sharp steps due to the two different Fermi energies in the two reservoirs (see 
Fig. I13|l . Such steps are measured by the Saclay group studying short metallic wires of type Fig.O) by attaching an 
extra tunneling diode on the side of the wire. (See Pothierl997 Gouaam2000). 

In the ballistic limit the distribution function is very similar only the factors ao(r) and 1 — Q!o(r) are replaced by 
the geometrical factors and 1 — determining the solid angles in which the ballistic electron arrives from the 
left or right reservoirs to the point r. 

The electron distribution function H126(l obtained for le\ does not result in any current in agreement that the 
resistivity is infinite in this limit. For finite 4i <SC d the factor ao(r) must be replaced by a momentum dependent one 
Q;p(r), and the current is due to da-pir) = a-p{v) — ao(r), which has a strong asymmetric momentum dependence and 
(5Q!p(r) ^ as le\ 0. 

The (5ap(r) will be determined by the very elegant theory of Kulik, Shekter and ShkorbatovfSJiiiJiiSSi who have 
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extended the treatment of the balhstic hmit to the diffusive one treating the scattering on phonons and defects 
which are large in space. The strong elastic scattering on defects is combined with a weak inelastic scattering. The 
limit considered is specified as (leiknY^^ ^ d where d is the size of the contact and ^oi and kn are the elastic and 
inelastic mean free paths of the electrons, respectively. Introducing the corresponding relaxation times, Td and r^n 
the inequality can be rewritten as (Tciri„)^/^ 3> That means that an electron can diffuse out from the contact 
region of size d with small probability of inelastic scattering, thus double inelastic scatterings can be neglected. In 
the following Tin is related to the relaxation time due to TLSs. 

Also in the diffusive limit the current correction due to the TLSs is derived by solving the Boltzman equation 
(|65|l with an appropriate collision term. In this case the collision integral /(p, r) has two contributions as the elastic 
impurity part /ci(p,r) and the contribution of individual TLSs at position R, /tls.rIPjT) thus 

/(p,r) =/ei(p,r) + ^/TLS,R„(p,r), (128) 

n 

where the sum is due to different TLSs, but in the following it is assumed that the concentration of the TLSs is so 
low that first a single TLS at position R is considered and in the final result the contributions of the different TLSs 
are additive. 

The kinetic equation can be arranged as 

vp^^ - eE^I^ Up, r) = /tls,r(p, r). (129) 
or op 

First the left hand side is treated in the diffusive limit and the right hand side is considered as a weak perturbation. 
The impurity part of the collision integral is 

''='(P'^)=(2i)^ / ^^^;'(/p'W-/pW)' (130) 

ep=ep/ 

where W^™? is the elastic transition probability, is the velocity at momentum p' perpendicular to the equienergetic 
surface for which the integral is performed. The electric field can be expressed by the electric potential $(r) as 
E = — Vr<I'(r), where $ is determined by the neutrality condition as in the ballistic case. 

The distribution function /p(r) satisfies the following boundary condition: At the boundary, E between the insulator 
and the metal we assume mirror reflection 

/p(r e E) - /p(r e E), (131) 

where an incoming electron with momentum p is reflected with the momentum p. Furthermore, very far from the 
contact |r| oo the equilibrium distribution is recovered with chemical potential /i, thus 

/p(|r| ^ oo) = /o(£p) = ^^^^ , (132) 

g fe^T 

and satisfies the charge neutrality. 

The distribution function in the absence of the collision term due to TLS is denoted by fp'' (r) and the first order 
correction due to the TLS is fp\r), thus 

/p(r) = /r« + 4'Hr) + ..., (133) 

where the higher order corrections are neglected. 

The kinetic equation for fp\r) contains the static impurity contributions, thus 



^p^^-^E(0)^^-Wr(r))=0, 



and /p^-* (r) satisfies the equation 



(134) 



(135) 
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J« = -2e / dS / ^^/W(r)vp = 0, 1, • • •), (136) 



where the electric field is also expanded as E = E^^^^ + E^^^ + • • •. The term /tls,r is linearized in the coUision, thus 
/p(r) is replaced by fp\Y). 

Similarly, the current is / = I^^^ + 1^^^ + • ■ •, where 

d^p 

(27rft)3 

and the integral with respect dS is taken on a dividing surface representing the point contact. 

According to the introductory remarks of this section, at finite elastic mean free path the distribution function can 
be expressed with a momentum dependent parameter, ap(r): 

/(")(r) = ap(r)/+ + (1 - a^{v))J^, (137) 

where the electrons arriving from far left (right) at the contact have distribution function (Jq), thus 

/o* =/o(ep-e$(r)±^). (138) 

As the collision term (I13U|I is linear in the distribution function the equation for /p*'''(r) is satisfied with ap{r) as well: 

Vp^-eE(o)^-/,(ap) = 0. (139) 

In the loi limit the solution of this equation is ao{r), given by (|127|l . At finite elastic mean free path ap{r) can 
be expanded as: 

ap(r) =ao(r) +5ap(r), (140) 

where Sap{r) is the first order term in the small parameter, The momentum dependent correction. Sap can be 
determined by a simple argumentation. At non-zero relaxation time the probability ap can be considered as the 
momentum independent probability ao taken at the position of the last elastic collision, that is: 

ap(r) = Q:o(r - TciVp). (141) 

The expansion of this formula in t^i gives: 

dap - -Tcivp^. (142) 

The value of ao is determined by the potential, which drops in the contact region, thus ^ ~ ^, and Jap - ^. 

The value of Sap can also be obtained by inserting H14(J|) into the Boltzmann equation for ap ni39|) . After neglecting 
the higher order terms in the small parameters ^ and |^ the following simple formula is achieved: 

Vp^-/ci(ap) = 0. (143) 



integral is expressed as: 



For isotropic scattering VFp"^'^ = ^^-p' holds, the relaxation time approximation is appropriate, thus the collision 



IJa )- f ^p^i-P - (144) 

which provides the relaxation time, t^i due to the elastic scattering by the impurities. Inserting this result into H143|) 
the formula (|142() for Sap is regained. 

The next step is to determine fp^ using the equation H143|l for ap. The boundary condition for ap(|r| — > oo) — 
Q{—z) as very left the electron comes from the left contact where the electrons are in equilibrium thus 

/W(|r|^oo) = 0. (145) 

Furthermore, the collision term contributes to the distribution function only if the system is out of equilibrium, thus 
/p^"* is proportional to the electric field. Therefore, the second term — eE*^"^^^- in equation H135|) can be dropped as 
it is proportional to |E^°)p. 
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The correction to the kinetic equation is determined by Eq. (|135|l where after hnearizing each term in the field the 
Ansatz with the function Xp(i") 



/W(r)=-e$W(r)^ + Xp(r) (146) 



f(0) 

can be used. The terms proportional to E^^^ ~ —d^^^^/dr cancel out in Eq. p35l) . Thus 



- ^ci(Xp) - /TLS,R(/r(r)), (147) 

with the boundary condition Xp(|i"l ^ "o) — ^'^d Xp(r G S) = Xp(r ^ S)- As it will be shown below the change in 
the distribution due to the collision term containing the interaction with the TLS can be expressed as linear response 
with the aid of a Green's function gpp>{r,r'): 

Xp(r) = J dh'J d3pW(r,r')/TLS,R(4?^(r')), (148) 



where 

dgpp'jr, r') 
P dr 

and 



+ /ci(gpp'(r, r')) = --^(P - P')^(i- - r'), (149) 



/ei(gppKr,r')) = J ^1^™?: {^p'V (r, r') - 5pp' (r, r')} , (150) 



with the boundary condition gpp'(|r| oo,r') — and 5pp'(r e S,r') = gpp'(r g 5],r'). As 5pp'(r, r') is closely 
related to the correlation function (np(r)np/ (r')) the time reversal symmetry can be applied as 

ffpp'(r,r')=g-p'-p(r',r). (151) 

The function Gp(r) is defined as 



Gp(r) - / dS'^^v^p,gpp,{r,r'), (152) 



where the integral with respect to dS' is a surface integral of variable r' taken on the orifice. This function satisfies 
the equation 

BC 

Vp^+/ci{Gp(r)} = -«;5(^), (153) 

which follows from the equation (|149|l defining gpp'(r,r'). The boundary conditions are the usual ones Gp(|r| — > 
oo) = and Gp(r g E) = Gp(r g S). This equation for Gp(r) is closely related to Eq. H143() for ap(r) introduced 
earlier, but without the inhomogeneous term. (The sign in front of the collision integral can be changed by replacing 
p — !• — p and taking the relation v_p = — Vp into account.) By comparing the two equations the expression 

Gp(r) = l-a_p(r) -e(z) (154) 

is obtained, which automatically satisfies the boundary condition for Gp(r). 

Using the expression H154() for Gp(r) the current can be calculated and only the p-dependent part is contributing, 
as the total number of electrons at position r is not changed by the collision. Therefore: 

yd3p/TLS,R(4°)(r))=0. (155) 

The kinetic equation 1)1471) for Xp can be checked using equations H148|l . (|149|l . (|15U|I and inserting twice the identity 
(|151|l . The next task is to calculate the correction to the current due to the inelastic scattering on the TLS. Using 

the expression p46|) for /p^^ (dropping the first term which is even in momentum p in leading order): 

I^'^ = JdsJ d^p v;xp{r), (156) 
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where the variable r is taken over the contact area, S. Inserting the expression of x the result can be rewritten as 

I^'^ = / d^d'P «P / d^P'd^r'5pp,(r,r')/TLS.R(4?^(r')). (157) 



That result can be further rewritten in terms of Gp(r) as 



(158) 



Using Eq. H154|l for Gp (r) the momentum independent terms do not contribute due to Eq. 1)155(1 . The remaining term 
contains 5ap which was determined earlier and 6a^p = —6ap holds. Furthermore, /TLS,R(Pjr) — 6(r — R)/tls,r(p) 
and in this way the final expression for the current is 



(159) 



where Sap is taken at the position of the TLS, R. Because of the odd parity of 6ap only the odd part of the collision 

integral contributes. Using the equation (|137|l for fp \r) the collision integral for the inelastic scattering can be 
written similarly to the ballistic case fEa. I71|l : 



m,TLS,in 



(P) 



{2nhf 



de'dfip/W^p/ 



"p'/o(£' + ^) + (l-«P')/o(£'-^) 



n-d{e' ~e-E)+ n+S{e' - e + E) 



"p(l - hie + ^)) + (1 - «p)(l - /o(e - ^)) 



P £ 



p' ^ p e' ^ e 



(160) 



is obtained, where in the second set of brackets the variables are changed as £ ^ e', p ^ p' and all the energy 
variables e and e' are shifted by — e<&(r) thus $(r) drops out. Now using Eq. ((14()|l ap = ao + 5ap and keeping only 
the odd part of the collision integral, after changing the variables the current is 



jd) 



2e 



An ^ in 



-Woo' 



6ap' 



eV eV 



"o(l - /o(e + ^ )) + (1 - ao)(l - hie - ^)) 



«o/o(e' + ^)) + (1 - «o)/o(£' - ^) 



eV eV ' 

h{s + —)-fo{s-—) 



[n_(5(e' -e-E)+ n+S{e' - e + E)] 
p' ^ p e' ^ e 



(161) 



where ao and dap are taken at the site of the TLS (r = R). 
In this expression two squared averaged matrix elements occur 



w 
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1 

(R) 



dQp 

An 



An 



^{Sapfwpp, = t\ 



dr^p 

An 



df7p/ 



[^f'^Sa,Sa-,=r!J^[^ 
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Wpp' 
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m dr I , 



2 "^F 
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dao 



dr 



dQp 

An 
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An 



p dao \ f p' dao 
m dr 



r=R 



dr 



Vpp' 



(162) 



r=R 



where the expression H142f) for Sap was used. Furthermore wi > \w2\ holds as Wpp' — Wp'p and 



p da 
m dr 



p' da 
m dr 
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p da 
m dr 



p' da 
m dr 
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Integrating over the variable e' the following expression is obtained 

/(I) ^ -^^/^r / d£ { foie - ^) - Me + 'Z)) X (163) 



1 + (n_ — n+) 



«o(/o(e + i?+^)-/o(£-i? + ^)j +(l-ao) (fo{e + E-^)-fo{e~E-^) 



where the geometrical factor Kfl is introduced for the diffusive limit as 

(R) (R) 

The expression ()163|l can be brought to a form which is comparable to the result H10n|l in the ballistic limit. The key 
of the algebra applied is in the integral 



de /o(e)(l - foie + a)) = J de foie + fe)(l - foie + a + b)) 



the variable can be shifted by an arbitrary energy b as the contribution comes from the neighborhood of the Fermi 
energy and the bandwidth cutoff does not play any role. The result of a lengthy algebra is 

/(I) = -y^^r J de foie) {{foie ~ E ~ eV) ^ foie -E + eV)) n_ + (/o(e + E - eV) - /o(e + E + eV)) n+} . 

(165) 

This expression coincides with the one obtained in the ballistic limit (|1UU|) after replacing the geometrical factor 
Kn, K^^. Therefore, the final results in the baUistic limit ((TUTl [TM [TUl^ are also valid in the diffusive limit, just 
the geometrical factor, Kn is different. 

The elastic contribution can be obtained by inserting E = Q and considering different scattering strengths for the 
two states of the TLS in the diffusive limit as well. Therefore, the ballistic result for the elastic case (|106ll is also valid 
in the diffusive regime with the replaced geometrical factor. 

The next step is to calculate n+. In the present case whether the electrons come from the left or right reservoir is 
determined by ao or (1 — ao), thus the previous results in the ballistic regime (|117l 11151 1118|l can be adopted after 
changing 

^ - ao. (166) 

The momentum dependent correction 5ap to ao contributes only if the momentum dependence of w-ppi is kept but 
that is small in case l^i <C d. 

Concluding, in the diffusive limit the shape of the current correction due to the scattering on a TLS is the same as 
in the ballistic case, the only essential difference is in the geometrical factors, which are independent of the strength 
of the interaction between the electrons and the TLSs. In the ballistic case i^R, is in the range of unity if the TLS 
is situated in the contact region. However, contains other quantities proportional to (tci^)^(^2.)^ ^ 'ci(^?')^ 

which indicates how much a is changed in space on the scale of the impurity scattering mean free path. We have seen 
that ^ scales with the inverse of the contact diameter, d, thus 

Kr-{^^)' (167) 

in case where the TLS is nearby the orifice. Far from the orifice (|r| — > oo) ^ — > 0. The amplitudes of the change in 
the conductance for the important TLS (being close to the contact) are scaled down by 

which can be a very small factor in the dirty limit. The largest contribution arises from the TLSs for which |^ is the 
largest, thus for the TLS just nearby the edge of the orifice. (See Eq. IS^ 

The change of the geometrical factor due to the impurity scattering length is similar to the case of electron-phonon 
interaction, where integration is taken according to location of the interactioni^^iiiiJiiSS The if-factors for electron- 
phonon scattering are given in (|67l I68|l for the ballistic and diffusive case, respectively. In the phonon case the 
amplitude of the signal in the diffusive limit is reduced only by (-j-)- 
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C. The case of time-dependent applied voltage 

Up to this point the scattering on TLSs in point contacts was investigated as a function of bias voltage, which 
provides information about the excitation spectrum of the TLS, and also about the nonequilibrium distribution of 
the TLS states. A further possibility to study the dynamical properties of the TLSs is the study of the response to a 
harmonic excitation. Let us consider that a voltage of 

V{t) ^Vo + Vi cos(wi) (169) 

is applied on the point contact, where the DC bias voltage, Vq is varied on the scale of the excitation energy of the 
TLS to cover the spectroscopic peak at e\V\ — E, and also the background tail at e\V\ > E. The amplitude of 
the harmonic excitation is considered to be much smaller than the excitation energy, eVi <C E. Furthermore, it is 
assumed that the energy fno is much smaller than the excitation energy of the TLS, thus the quantum nature of the 
radiation should not be taken into account. In this case the current flowing through the contact can be expanded in 
second order in terms of the time dependent voltage contribution: 



J(Fo + V, cosM) = ^(^") + i 



V-f H 

-V dV 



FrCOs(c.t) + -^ 

-Vn 



V^cos{2ujt)+..., (170) 

V=Va 



where the equation cos^(</5) = (cos(2(p+ l))/2 is used. The voltage dependence of the harmonic response, I'^ provides 
the differential resistance curve, dl/dV. The second harmonic response, P'^ directly provides the point contact 
spectrum, i.e. the second derivative of the liV) curve. The DC component, which can be measured as a time average 
of the response signal: /° = (/(Vq + Vi cos{u)t))) also shows the second derivative of the I{V) curve after subtracting 
the current without harmonic excitation, /(Vq). Further on this term is referred to as the DC shift signal, denoted 
by 5P. 

The above considerations are, however, only valid if the system is fast enough to follow the time dependent voltage. 
In the case of two level systems, the dynamical behavior is determined by the kinetic equation (|115|l : 

^=A-n+-{A + B), (171) 

where the coefRcients A and B are given by Eq. Hllt)|) . Both of the coefficients depend on the actual voltage V{t) at 
given t. For simplicity the notation A + B — B \s introduced, which is the inverse of the relaxation rate of the TLS 
{B = r^^). For slow two level systems the value of the relaxation time can be estimated using Eq. (|115I91I70|I as 
follows: 

Bc,—{2H'{QoVn' (172) 

Inserting ^ = = 1/2, = 1 meV, and goV ~ 0. lYi'^'^'^''i???'^-^i=^'^'^''i???c into this equation a relaxation time of 
T ^ 10"^" s is obtained. 

If the alternating voltage is faster than the relaxation of the TLS {lut ^ 1) then a more complicated procedure 
is necessary to determine the time dependent response of the system. This phenomenon, however can be used to 
trace the typical relaxation time of the TLS by studying the response signal as a function of the frequency uj. The 
first calculation with time dependent applied voltage was performed by Kozub and Kulikii^saii^i22& In the next part 
a general solution is given for the time dependent occupation number n^{t) at arbitrary frequency and amplitude of 
the harmonic excitation signal. After that various limits are treated, where the formulas are essentially simplified. 



1. General solution for n+{t) 
The kinetic equation l|171|l is a linear inhomogcneous differential equation which has a general solution 

t _ / * t' _ N 

f df' B{v{t')) / r r dt" B{v{t")) , 

n+{t) - e ■\Ci+ dt' A{V{t'))J^ , (173) 



where the coefficient Ci is determined by the initial condition at time t = Q. According to the definitions given by 
Eq. H116|) A>Q and _B > hold, thus the exponent in the last term of Eq. (|173|l monotonically increases with time. 
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This means that Ci is essential only for the transient behavior and for long time it can be dropped. For periodic 
voltage A and B are also a periodic function of time: A{V{t)) = A{t) = A{t + T^), B{V{t)) = B{t) = B{t + T^). The 
time integrals can be divided according to the periods as: 

Jdt... = J dt+ j dt+ ...+ j di+y^ di, (174) 

T„ (ti-1)T„ nT^ 

where t = nT^ + f, n is an integer and < f < T^^. In this way: 

t 

j dtB = nAQ + Q{t), (175) 



where 

AQ = j dtB (176) 



and 

i 

Q{t) = j dt' B{t'). (177) 



Using these identities, one gets: 

* - „ „ n-l * 

/ dt' A{t')J> = e"'"^"^ / di' A(i')e^(*'' + e"^« j dt! A{t')e'^^''\ (178) 



"'=0 



The final result is: 



= ^ ( ^"i^j^t' A(0e'3(*') + Jdtf A(t')e<^^^'^ \ , (179) 





which gives the following stationary solution in the long time {n^ 1) limit: 



(180) 



Turning to dimcnsionlcss variables like (p — 27ri/T[^, and (p = 2'Kt/Ti^ and also new notations, a{ip) = A(V{t)) /2Tr, 
b{ip) = B{V{t))/2TT, q{ip) = J b{ip')dip' = Q{t)/T^ and Aq = g(27r) = AQ/T^ one gets: 



^tJ,_-^ T^ j V a{<fi')e^-''(^'^ +T^j d^' a{v')e^-'^^'^') . (181) 
/ 

It is important to note that q{(p) is proportional to the relaxation rate 1/r and a{(p) to the excitation rate of the 
upper level, which are independent of the frequency uj. 

In the following two simple limits are treated: wr — > {T^/t — > oo); and — > oo {T^^/t — > 0). 
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(a) (b) (c) 

FIG. 14: Three different cases for the harmonic excitation with respect to the excitation energy of the TLS. 



a. Low frequency limit (lot <C 1 J Assuming that Tt^q{(p) 3> 1 {Tuj/t 3> 1) the first part of the right hand side 
of Eq. H181|l disappears as q{f') < Ag = <z(27r). (Note that q{ip) is a monotonous function.) The contribution to 
the second part also disappears as far as q{(p') — q{<p) < is finite. In the vicinity of the upper fimit of the integral 



[if = if), however, q{ip) can be expanded as q{(p') — q{ip) — 



V). where 



b{ip'), thus 



T,, 



(182) 



where the term arising from the lower limit is dropped as it is exponcntionally decreasing with T^^ 
result for the limit — )■ cxd is: 



n+{t) 



a(^) _ a{^) _ A{V{t)) 
b{^) b{^) B{V{t)) 



oo. The final 



(183) 



That result also holds for ^ 1, where the condition T^q{^) 1 is not satisfied any more. Without making a 
rigorous proof, it can be argued, that <C 1 is at the beginning of the time period (p = Q. In the long time limit 
(t 3> TL) the time can be shifted , therefore, that time is not special, thus the result (|183|l must be generally valid. 

In this way the static result is recovered with A and B taken with the actual voltage at time t. Therefore, at 
arbitrary time the system is in the stationary state, and the the expansion of the current in (|170|l is valid. 

b. High frequency limit (ujt 3> 1 j Taking only the first part of the right hand side of Eq. (|181f) contributes. 

The denominator can be expanded as e^"'^'' — 1 ~ T^Aq and the result is 



71+ (t) = 



27r 





^ M 

I dtp b{ip) 





(184) 



where (...) stands for an average for a time period and the definition of Aq is used. As the TLS reacts slowly on that 
time scale, thus it reacts on the average value of the parameters A and B. Note that {A)/{B) ^ (A/B). 



2. Dynamical conductance 

In order to calculate the time dependent current first the occupation number n^{t) must be calculated, which is 
not an easy task in the general case as the formula given by Eq. (|181l) contains double integrals. The next step is to 
insert the obtained n+{t) into the equation for the current (|101|l . If one is interested in the average current (DC) and 
the harmonics of the current, it is not enough just to determine the the averaged and its harmonics. 

That can be demonstrated easily at T = by Eq. H1U3|) . At e\V\ > E the expression contains a term, which is 
independent of the occupation number. This causes a nonlinearity in the I(V) curve in every case, regardless of 
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the voltage dependence of the occupation number. For instance, if the TLS was in thermal equilibrium with the 
bath then n+ = would hold at any bias voltage. In this case the I{V) curve is linear both at e\V\ < E and at 
e\V\ > E, but it has a breaking point at e\V\ = E. If the alternating voltage is smaller than E in the whole time 
period (e|Vo + ViCos{Lut)\ < E), then the response to the excitation is completely linear, and no higher harmonics 
are generated. This situation is demonstrated in Fig. 114b. Similarly, the response is linear if e\Vo + V\ cos(ijji)| > E 
at any time. However, if the condition e|Vb + Vi cos(cL;i)| < -B is only valid in a part of the time period, the system 
gives a rectified response to the harmonic excitation, thus higher harmonics are also generated. This case is shown in 
Fig. 114b . The TLS can be excited only in the part of the time that is demonstrated by the shadowed areas, where 
the conductance is reduced due to the possibility of backscattering. In Fig. [Tlh a third case is presented, where the 
amplitude of the alternating voltage is so large, that rectification occurs both at eV — E and at eV = —E. In this case 
the contributions from the both polarities work against each other, and the higher harmonic generation is reduced. 

In a general case both the nonlinearity of the occupation number n^{V) and the above discussed intrinsic non- 
linearity of the I{V) curve cause higher harmonic generation. However, if the period time of the excitation is much 
shorter than the characteristic relaxation time of the TLS then the occupation number cannot follow the alternating 
voltage, thus it sets at a constant steady state value. In this situation the higher harmonic generation is still present 
due to the intrinsic nonlinearity of the I{V) curve. 

In the following section these features are discussed in two special cases. In Sec. IIV C 31 the T = the limit is 
treated in the high frequency limit {lot — > cxd), where the rectified current and the harmonics can be easily calculated 
at arbitrary value of the amplitude of the excitation, Vi. These results are also valid at finite temperature if eVi S> kT. 
In the second case (Sec. IIV C 4|) an expansion with respect to Vi is applied. This approach is valid if the amplitude 
of the excitation is small {eVi < kT), however in this limit the frequency dependence can be investigated in a wide 
range of frequencies. 



3. Rectification at zero temperature in the lot oo limit 

The rectification due to a TLS in the contact can be nicely demonstrated at T = in the lot oo limit. We 
have seen that at lot — > oo the occupation number becomes time independent which makes it possible to perform 
the calculation analytically. According to Eq. 11841 lim n+ = , where (A) and (B) are the time average of these 

coefficients. Inserting zero temperature into the formula (|116|l the voltage dependence of A and B can be written as: 

if e\Vo + Vi cos{Lot)\ < E 

K{e\Vo + Vi cos(wt)| - E) if e|Vb + Vi cos(wt)| > E 



m - j-iE + 2A{t)) , (185) 

Therefore, in case of amplitudes for which e|Vb + Vi cos(w<)| < E holds for the total period of time the occupation 
number is zero. If that is not satisfied, then the actual voltage at a certain part of the time can excite the TLS out of 
the groundstate. Two cases must be distinguished: (i) where Vb + Vi > E but Vo — Vi> ~E; (n) where Vq + Vi > E 
and Vo-Vi < -E. 

Figure 1151 show the actual voltage as a function of time using ip = Lot in these two cases. The shadowed areas 
indicate the regions where e|Vb + Vi cos{ujt)\ > E. The phases, where the shadowed areas are ending are denoted by 
ipo (in case (i)) and (p = tt — ipQ (in case (ii)). All of the following results for case (i) can be applied also for case (ii) 
by inserting ipg = 0. These phases can be expressed by Vb, Vi and E as follows: 

E~eVo ,_, , -E-eVo 

cos(v3o) = 77 — ; cos(.^ = - cos = . (186) 

eVi eVi 

The time average of A is calculated as: 

- f ? } 1 

{A) = -l / [e{Vo + Vi cosv?) ~E]d<f+ / [-e{Vo + Vi cos(p) - E]d<f\ = 

I ifio ) 

= -e [l^,(¥'o - <^^) + Vi (sin ipo + sin ip*)] - -E{ipa + </?*), (187) 

TT TT 
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case (i) case (ii) 

FIG. 15: Illustrations for the calculation of rectification. 



and the occupation number is written as 



(A) 



After having the result for the inelastic current correction can be calculated using Ea. I1U3I At the points, where 
e|Vo + Vi cos(a;t)| = E the time derivative of the current is discontinuous, which results in higher harmonic generation. 
The second harmonic generation with cos(2wt) and sm{2ujt) and also the shift of the DC signal can be easily calculated: 

— / I{V{ip))dip = 2eeoW-f^fl — (1 - 2n+){smipo - (pocos(po - sinipQ + ipqCosipq) (189) 

ZTT J TT 



J" = 





27r 



jZs 



1 f 

— / I{V{(p))cos{2(p)dip = 2eglwKR^{l~2n+){sm^^o-sm^^o) 



^ J I{V{^))sm{2^)d^ = Q 



(190) 



(191) 



The results for the cos(2a;i) are plotted for different values of the excitation amplitude, Vi. In this case the width 
of the peak at e|Vb| = is determined by the amplitude Vi, and these results are also valid at finite temperature if 
kT<^eVi. 



4- Expansion with respect to Vi 

At finite frequencies the occupation number is already time dependent (n+(i)), thus the above considerations 
cannot be applied. A calculation for intermediate frequencies can be performed in the eVi <JC kT limit. In this case 
the integrals in the general solution for n+ fEa. I173|l are calculated after expanding the coefficients A{Vq + V\ cosujt) 
and B{Vq + V\ cosujt) with respect to Vi. The expansion for A is: 

A{V{t)) = Ao+ A'qVi cosujt + cos^ ojt + ... (192) 

and similarly for B. The derivatives are taken at y = Vq. In the general solution for first the quantity 

must be calculated. The non-oscillating part of the exponent is linearly increasing with time, thus the expansion can 

be carried out only for the oscillating part assuming that BqVi/uj <^ 1 and Bg Vi /Sui ^ 1: 

Jdt'B(t') / , sinoji 1, /,9 o/sinojA^ 1 -r, „ sin 2a; A 

^e^^-il + B,V,-- + -iB,fV,^(--\ +-,BoV,^^], (193) 
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FIG. 16: Second harmonic generation at zero temperature in the cut 
of the excitation signal. 



oo limit. The curves are plotted for different amplitudes 



where 



(194) 



The shift of the relaxation rate by a term proportional to Vi is due to the nonlinear behavior. The next step is 
to calculate the integrals / dt A{t)e'^'^*\ All the terms obtained are proportional to the exponential, e^*, which is 
cancelled out by the factor e~'^* in the expression H173|l . In the remaining part the higher order terms 0{Vi) are 
dropped. The next tedious task is to collect all the terms in the coefficients of the expansion. The occupation number 
can be written in the following form: 



n+ = n"f_ + Vi{n\_ sinujt + n\ cosujt) + V^{5n°_^ + n'^_^ sm2Lut + n^_^ cos2ujt) + 0{Vl), 



(195) 



where n^J. is the stationary occupation number &t Vi = {nP^ — Aq/Bq), and define the amplitudes of the 
harmonic response to the excitation, (5n^|_ is the shift of the stationary value due to the nonlinear behavior, while n]^ 
and n"^ define the amplitudes of the second harmonic generation. The values of these coefficients are given by the 
following equations: 



5nl 



,2s 



- AoB'o 



A. 

4/3 
iA'aUj^ 



LU/3 



Ao(So)2 A„Bou;^ + AoiBo)^P 



Af3Bo 2^2/5 2u;2(/32 + 

AAoiBo)^ + AoB^p + 4A'oBal3 



^0^0/3-^0(^0)' AoB„ 



8tj(/32 - 
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2w(/32 + uj) 



^0^0^' ■ 



8cj/3 

^(5;)2/3 Ao(Bo)2 
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2tj2(/32 
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(196) 
(197) 
(198) 
(199) 
(200) 



The conditions for the expansion have excluded the narrow frequency region around lu = 0. Accordingly, the above 
results contain singular terms in uj, which diverges as w — > 0. On the other hand, it has been seen that the a; — > 
limit is well defined for 71+, and indeed, it can be shown that the singular terms cancel out in the low frequency limit. 
The expression for the current correction can be expressed similarly to l|195|) : 



1 = 1" 



Vi {r sin Lut + I" cos Lut) + {61° 



sin 2ujt + I^" cos 2ujt) + 0{V^), 



(201) 
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FIG. 17: The left panels show the voltage dependence of DC current shift (Slf„) and the second harmonic terms (ifn, if^) 
as a function of bias voltage at different excitation frequencies in the inelastic case. The right panels presents the frequency 
dependence of the corresponding terms at a fixed voltage of eV — E. The temperature is kT = O.IE. 



in the inelastic case the current correction given by Eq. (|1()1() can be written as 

Iin=n+D + n^C ^n+D + C, (202) 

where D = D — C. These must be expanded similarly to equation H192|) given for A. After this step the coefficients 
in H201() can be given as follows: 

I-n = KDo (203) 

/f„ = n%Do + nlD'o + C'o (204) 

51°, = 6nlDo + + + ^n^D, (205) 

= nlWo + ^n^+D, (206) 

= nl^Do + \C'^ + + ^<^o (207) 

These results can be analytically calculated, and for demonstration the experimentally important quantities {Sl'^, 
I^'^, and /^^) are plotted in Fig. El The left panels show the voltage dependence at different uj values, while the 
right panels demonstrate the to dependence at the fixed bias voltage value of eV — E. At low frequencies, the TLS 
occupation number can follow the alternating voltage, thus both P'^ and 51° are proportional to the second derivative 
of the I{y) curve (see Ea. ll7()|l . In this limit there is no phase shift, and thus = 0. In accordance with the static 
results, the curves show a negative peak at eV = E, and a positive tail at higher voltages. In the high frequency 
limit the TLS cannot follow the excitation, thus the occupation number sets at a constant value. Due to the intrinsic 
nonlinearity of the I{V) curve, even in the high frequency limit a pronounced negative peak is observed at eV — E, 
but the positive background disappears. (At zero temperature the amplitude of the negative dirac delta peak at 
eV = E in P'^ is exactly 2/3 times smaller in the a; oo limit than in the lo ^ limit. This factor of 2/3 is well 
demonstrated by the w dependence of P'^{eV = E) in Fig. El as well.) Around the characteristic frequency of the 
TLS {lot ~ 1) a significant phase shift is observed, as demonstrated by the curves for . 
In the elastic case the different current components can be expressed by using Eq. IjlOTf) : 

II, = -(^+-^-)Fo< (208) 
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FIG. 18: The left panels show the voltage dependence of DC current shift (Sl'^i) and the second harmonic terras (I^i , lli ) as a 
function of bias voltage at different excitation frequencies in the elastic case. The right panels presents the frequency dependence 
of the corresponding terms at a fixed voltage of eV = E. The temperature is kT — O.IE. 
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In Fig. ^] the second order terms are plotted for the elastic case. The nonlinearity of the I{V) curve of the elastic 
contribution is basically determined by the nonlinearity of n^{V). Above the characteristic frequency of the TLS the 
occupation number cannot follow the excitation, thus both /^j" and I^f are suppressed. The shift of the DC component 
((5/"j) shows the spectroscopic peak at eV = E even at lot oo, just the amplitude is reduced by ~ 30%. Around 
the characteristic frequency of the TLS a significant I'^y contribution is observed. 

The dynamical properties of the system are also demonstrated by the phaseshift in the second harmonic contribu- 
tions: I^f cos{2u;t) + I^f cos{2u;t) — I^^^ cos(2cjt + (p). In Fig.^lthe phaseshift ip is plotted at eV ^ E both for the 
inelastic and elastic case. In the elastic contribution the phaseshift grows from to 90° between the ujt and the 
LOT — > OO limit. Contrary, in the inelastic case the phaseshift is zero in both limits, and even at lot ~ 1 a only minor 
phaseshift of ~ 10% is observed, as the nonlinearity is basically related to the intrinsic nonlinearity of I{V) curve and 
the nonlinearity of n_|_(l^) plays a less important role. 



D. Conclusions 



In this section the theory of the scattering on slow TLSs in point contacts was reviewed. In Sec. IIII Cl it was seen 
that the electron distribution functions and the electric potential around a ballistic contact can be determined by 
solving the homogeneous Boltzmann equation for a point contact geometry. The solution yields a voltage independent 
conductance, which is known as the "Sharvin conductance" fEa. lt)3|) . The effect of the scattering on a slow TLS near 
the contact can be considered by perturbation in the collision term of the homogeneous Boltzmann equation. The 
electron-TLS interaction may result in backscattering processes, where an electron that has already crossed the contact 
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FIG. 19; Phaseshift in the second harmonic current correction at eV = E for the elastic and the inelastic contribution, 
respectively. The temperature is kT — O.IE. 



is scattered back through it. These processes give a strongly voltage dependent correction to the current, for instance 
at zero temperature backscattering can only occur if the voltage bias exceeds the excitation energy of the TLS. 

The nonlinear current correction due to the inelastic scattering of the electrons on a single TLS can be explicitly 
given, if the occupation number of the upper level n_|_, the electron-TLS coupling strength w, and the geometrical 
factor i^R, are known. It has crucial importance that the interaction of a TLS with the neighborhood is dominated 
by the electron-TLS scattering, and the phonon contribution is significantly smaller because of the relatively small 
phonon phase space at energies much smaller than the Debye energy. Therefore, the excited TLS mainly relaxes its 
energy through the electrons, which have a nonequilibrium distribution near the contact. As a consequence, the TLS 
occupation number is also not the equilibrium one, and it is strongly voltage dependent. The value of n+{V) can be 
determined by a simple dynamical equation using Fermi's golden rule. Due to voltage dependence of n_(_ a nonlinear 
current correction arises from the elastic scattering of the electrons on the TLS as well, provided that the scattering 
cross sections are different for the two states of the TLS. Both for the elastic and inelastic scattering the point contact 
spectrum (i.e. the second derivative of the I —V curve) shows a sharp peak at the excitation energy of the TLS, and 
a background tail at higher voltages due to the nonequilibrium occupations. 

The current correction due to the scattering on a single TLS can also be determined for a diffusive point contact. 
In this case the equations are more difficult, as the collision integral due to the elastic scattering on defects should 
also be included in the Boltzmann equation. Surprisingly, the lengthy calculations give exactly the same results as in 
ballistic contacts, just the geometrical prefactors are different. 

The response to a time dependent bias voltage was also reviewed. Due to the nonlinearity of the I — V curve the 
application of a small AC modulation beside the DC bias causes two important features in the response signal: (i) 
a second harmonic signal is generated, (ii) the average value of the current is shifted, which is referred to as the DC 
shift signal. If the frequency of the AC modulation is much smaller than the relaxation frequency of the TLS and its 
amplitude is small compared to the energy splitting, then both of these signals are proportional to the second derivative 
of the I — V curve. On the other hand, if the modulating signal is fast compared to the relaxation of the TLS then the 
occupation number n+ cannot follow the variation of the bias voltage, which has important consequences. Previously 
it was believed that in the high frequency limit both the elastic and inelastic contributions from the scattering on 
the TLS are completely suppressed in the response signal. Our more detailed calculations have shown that the DC 
shift signal contains a frequency independent contribution due to the rectification at voltages in the range of the 
TLS energies thus it is not decisive for the relaxation dynamics of the TLS. We have demonstrated that both the 
elastic and inelastic contributions are almost independent of the irradiation frequency in the DC shift signal. A more 
appropriate quantity for these studies is the second harmonic generation signal. In this case the elastic contribution 
is completely suppressed and the inelastic contribution shrinks by a factor of 2/3 if the irradiation frequency exceeds 
the relaxation frequency of the TLS. We have also investigated the phase shift in the second harmonic signal which 
might be examined in subsequent experimental studies. 

There is another suggestion by Kozub and RudirJ^2£li^i22i to explain zero bias anomalies based on the adiabatic 
renormalization of the splitting of the TLS by electrons. That gives, however, a weak temperature dependence 
because the main contribution comes from the electrons far from the Fermi energy, thus that is not related to infrared 
singularities. The interaction between the TLS due to the Friedel oscillation was also considered and magnetic field 
dependence is suggested. However, it is hard to compare the theory with experiment in its present form. 

It is also interesting to note that TLSs and magnetic impurities interacting with metallic electrons show several 
similar features. Both could lead to logarithmic corrections known as the Kondo effect. The main similarity is that 
the spin S ~ 1/2 has also two states. Without external magnetic field or internal magnetic field characteristic to 
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FIG. 20: Resistance vs. time in copper nanobridges for T < 150K showing several types of behavior. The fluctuations range 
from 0.005% to 0.2% of the total resistance. The time scales depend on the temperature at which the fluctuation is observed, 
(a) A single TLF. (b) T wo independent TLFs. (c) Amplitude modulation, (d) Frequency modulation of one TLS by another. 
[Taken from lRaUslQSl . with kind permission of the American Institute of Physics] 



the spin glass systems the two states — ±1/2 are degenerate, just like a TLS without sphtting (A — 0). In that 
case the occupations of the two spin states ni/2 — n^xji = 1/2 and that cannot be affected by the apphed voltage 
on the point contact because of spin degeneracy. In Born approximation that gives a single peak at zero bias in the 
point contact spectra and no background occurs. The logarithmic Kondo corrections, however, broaden the 

peak which looks very similar to background in case of phonons and TLS, but that is not due to a nonequilibrium 
distribution£ia£iJa!i£li2Jia^ The situation is different, where the spin degeneracy is lifted by the magnetic field. In 
that case ni^2 ^-i/2j thus the occupation numbers depend on the applied voltage and background occurs even 
without Kondo corrections. 

In the following section the former experimental results are discussed and confronted with the model of slow two 
level systems. 



V. DISCUSSION OF THE EXPERIMENTAL OBSERVATIONS 



This review is aimed at collecting and summarizing the main results concerning the theory of the scattering on 
slow two-level systems in point contacts, and to compare those with the experiments. As a point contact probes the 
restricted area around the contact region, it can be used to investigate a few or even a single TLS located in the 
contact area. Therefore, some measurements can be directly used to study the nature of TLSs in the contact region, 
and to check the relevance of the model presented. 

In disordered metallic point contacts TLSs with a wide variety of relaxation times can be present. A spectacular 
sign of the presence of TLSs in a point contact is the telegraph fluctuation of resistance, caused by defects fluctuating 
between metastable configurations on an experimentally resolvable timescaleiSailSiSS Examples for such so-called two 
level fluctuators are shown in Fig. [201 Panel (a) shows the simplest telegraph fluctuation: the resistance is switching 
between two values, as discussed by considering the elastic processes (see Sec. lIVir2|) . while the switching itself is an 
inelastic process. In panel (b) two independent telegraph fluctuations are superimposed upon each other, while on 
panel (c) and (d) one two-level fluctuator modulates the amplitude or the average switching time of the other. The 
screening of these two-level systems by electrons causes a slow-down of the TLS motion which can be described by 
the renormalization of the tunneling amplitudeiM£ig£Sii£iia£ti^®'^''^'°'''^^"°i^9i This phenomenon can be used to estimate 
the electron-TLS coupling constant, as it is discussed in [Goldingl992. and .Zimmermanl99ll 

This review concentrates on slow TLSs where the relaxation time is much larger than the electron-TLS interaction 
time, but still the TLS is fast enough not to be able to resolve it as a telegraph fluctuation of the resistance. It 
was seen, that such "slow" TLSs can be studied by point-contact spectroscopy, a singularity in the PC spectrum 
(d'^I/dV^) is expected if the voltage bias coincides with the energy splitting of the TLS close to the contact. Figures 
[TUl and ^] demonstrated that the contributions of both the elastic and inelastic scattering exhibit this anomaly at 
eV — E. In the inelastic case this anomaly always corresponds to a maximum of the differential conductance, dl/dV, 
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which is reflected by a negative peak in the point contact spectrum d"^! jdV"^ at eV = E. In the elastic case a similar 
peak is observed, however its sign can either be negative or positive, depending on the sign of 7+ — 7" (see Ea. llU7|l . 

A similar peak is expected in the PC spectrum for fast Kondo-like TLSs, for which the electron assisted transition 
processes are not negligible (see Sec. lIIB'2|) . For the fast TLSs the anomaly appears as a minimum in the differential 
conductance around zero bias, and the width of the zero bias minimum is determined by the Kondo temperature 
Tk- This feature is shown as a positive peak in the d'^I/dV'^ curve at eV = fcrj^i?"^^???!^'^'?''^??^!^'"'?''^??^ That peak 
shows very strong similarity to the Kondo peak in the presence of magnetic impurities^SiSSSaiSSi but no magnetic 
field dependence is observed, and the peak is changing by annealing. 

The sign of the anomaly has crucial importance in distinguishing the slow TLSs from fast TLSs. Along this 
discussion the anomaly exhibiting a conductance maximum at zero bias is called "positive" zero bias anomaly (ZBA), 
while the conductance minimum is called "negative" ZBA. We note that the labelling of the sign of the anomaly can be 
confusing, as in some cases the voltage dependence of either the differential conductance or the differential resistance 
is plotted, in other studies the point contact spectrum is presented. Furthermore, concerning the PC spectrum, in 
theory the d"^! jdV^ curves have direct physical meaning (showing the excitation spectrum), whereas in experimental 
studies usually the second derivative of the inverse curve, d^I/dV'^ is plotted. Regardless of the way how the PC 
spectrum is plotted, the sign of the anomaly can be determined by comparing it with the sign of the contribution 
of the phonons. For the inelastic scattering on TLSs the increase of the bias voltages causes an enhancement of the 
back-scattering, thus this anomaly has the same sign as the phonon part of the spectrum. In contrast, for Kondo-like 
anomalies a resonant scattering is observed, which is restricted to the energy range of the Kondo temperature. In 
this case the back-scattering is enhanced at zero bias, and it decreases as the bias is elevated above eV ~ fcTjf , thus 
the anomaly has negative sign compared to the phonon contribution. In terms of this labelling, the fast TLSs should 
always show a negative anomaly, the inelastic contribution from slow TLSs shows up as a positive anomaly, while the 
elastic contribution has an indefinite sign. 

In the experiments different zero bias anomalies were observed. A basic distinction should be made between the 
experiments performed on pure metals, and those, performed on disordered systems, like metallic glasses. In the 
first case the temperature dependence of the bulk resistance does not show any anomalous behavior, but due to the 
magnifying effect of point contacts, the scattering on a few defects near the contact might be seen in the point contact 
spectrum. In this case the appearance of the anomaly depends on the way the contact is created, and usually the 
anomaly can be destroyed by annealing. In metallic glasses, however, a high degree of disorder is initially present, 
and the anomalous feature seen in t he PC spe ctrum is also reflected by the anomalous low-temperature behavior of 
the bulk resistivitjiSssisiSSi (see also lCoxl999|) . 

Several metallic contacts show anomalies with different signs depending on the way the contact is prepared. In the 
following we concentrate on a few typical and most interesting cases. 

A. Measurements on Cu point contacts 

The zero bias anomalies were widely studied in copper nanocontacts created with different methods. The first 
experiments by touching two Cu electrodea^JSiSiSS^SiS^^ reported a pronounced positive anomaly at ~ 1 meV [Fig. 121b]. 
The ZBA has a highly asymmetric shape, which can be well fitted with the theoretical prediction for the elastic 
scattering on a single slow TLS. The inelastic scattering could be excluded, as in this case the anomaly should show a 
basically symmetric peak with a background of opposite sign [see Fig.^]. Furthermore, as the ZBA has consequently 
shown a positive sign, a Kondo-like explanation could also be excluded. In some occasions further low-energy peaks 
were also observed [Fig. 121b]. which were argued to be due to the higher excitations of the TLS. 

Similar ZBAs were seen in the experiments of Keijsers et aliJ^SlS£Si22& on copper thin film mechanically controllable 
break junctions [Fig. 122b ] . These samples were created by vapor deposition, and judging from the residual resistance 
ratio and the intensity of the phonon peaks in the PC spectrum they were of low crystalline quality implying that the 
contacts were in the diffusive regime. However, the authors claim that the ZBA is not due to the intrinsic disorder 
of the sample, but it correlates with the amount of bending and stretching of the film. In these thin film samples the 
ZBA has consequently shown a positive sign, and it could be well fitted with the elastic scattering on a slow TLS, 
similarly to the above discussed measurements of Akimentko et al. 

The same work of Keijsers et al.^'^"^'^''^^^^^ has also presented measurements on single crystal Cu break junctions 
[Fig. 122b ]. In this case the well resolved electron-phonon interaction spectrum implies a large elastic mean free path. 
In these measurements the ZBA has consequently shown a negative sign. Here, the elastic scattering on a slow TLS 
is still a reasonable explanation, however, the Kondo-like behavior also has to be considered. 

Ralph et aliSaiEiii22^ have investigated nanofabricated Cu contacts. In this case the size of the contact naturally 
could not be varied during the measurements, but as a great advantage, the contact was "fixed", and thus a detailed 
study of the temperature and magnetic field dependence could be performed. In these experiments a well defined 
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FIG. 21: (a) Point contact spectrum (curve 1) and spectral density of noise (curve 2) for a copper point contact. The circles 
show the theoretical prediction for the elastic scattering on a single slow TLS. The inset shows theoretical curves for the point 
contact spectra and the spectral density of noise, (b) Point contact spectra for copper contacts showing several low energy 
peaks (marked by arr ows). The inset shows the possible transitions in a double-well potential leading to the peak structure 
observed. [Taken from lAkimenkol993l with kind permission of Elsevier Science.] 




FIG. 22: (a) Point-contact spectra of copper thin-film MCB junctions with resistances 1-1. 2f2, 2-220, 3-6457 and 4-40057 
(T = 1.3K). The curves have been shifted for clarity, (b) Point-contact spectrum for a Cu single crystal MCB junction 
{R = 30Q, T = 2K). The peaks at 17-18mV and 30mV are due to the electron-phonon interaction; the feature at bias voltages 
between and 5mV [positive for (a), negative for(b)] is attributed to interactions between electrons and TLS. Similar anomalies 
were also observed in Au, Cu, Ag, and Pt contacts. [Taken from lKeiiscrsl995i with kind permission of the American Institute 
of Physics.] 



negative anomaly was observed if the sample was quenched [Fig. I23j . and the anomaly disappeared after annealing. 
The detailed theoretical analysis of the curves (including scaling arguments)i!2ii2£i£*i^^* has shown very good agreement 
with the theory of fast Kondo-like TLSs, even if their existence is still an open question,^i£iiiaSaaj^ii^l£iii£'-20oib.Borda2^ 
The relatively high Kondo temperature {Tk ~ 10 K) can be obtained in the framework of recent theorji22l^s2222i only 
if the center is light (e.g. H or He) as the energy scale is the kinetic energyi^2£SitesJsi2224 

In all of the above experiments the ZBAs were attributed to a single or a few TLSs. In the measurements, where 
consequently a positive ZBA was observed both the inelastic scattering on a slow TLS and the Kondo-like behavior 
can be ruled out due to shape and the sign of the anomaly, respectively. Therefore, the clastic scattering on a slow 
TLS remains as an explanation, which indeed gives good fit for the curves. However, it should be noted that for 
elastic scattering the sign of the anomaly is expected to be contact dependent. A consequent positive sign might be 
explained by the following argument: the atom in the higher level of the TLS can be considered to be more far away 
from the ordered arrangement of the atoms, thus its scattering cross section is larger than that of the low energy 
state, which results in an enhanced backscattering as t he bias is incr eased. 

It is interesting to note that the positive ZBAs in lKeiisersl995l exhibit a pronounced contact size dependence: 
as the diameter of the junction is decreased the ZBA shifts to higher energy values [Fig. 122b]. The ZBAs showing 
Kondo effect in magnetically doped contacts have exhibited a similar size effect .^iSgSaJS The latter phenomenon was 
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FIG. 23: (a) Point contact spectrum for an unannealed Cu nanoconstriction. The negative zero bias anomaly at ~ 2 mV is 
attributed to the scat tering on fa st two level systems, while the structures at higher voltages originate from the scattering on 
phonons. [Taken from lR,alDhl994l with kind permission of the American Institute of Physics.] (b) Differential conductance of a 
Cu nanoconstriction at T = 100 mK. The curve measured in 6T magnetic field is shifted down by 20e^//i for clarity, and only 
a slight dependence on the field is observed, which supports the nonmagnetic origin of the anomaly. For comparison the inset 
shows the conduct ance of Cu constriction with 200 ppm Mn at 100 mK, showing a pronounced Zeeman splitting in magnetic 
field. [Taken from lRalphl992l with kind permission of the American Institute of Physics.] Similar zero bias anomaly was also 
observed in Ti nanojunctionsiiifiaaaaJiiffi 

explained by Zarand and Udvardi as a result of density of states oscillations nearby the surface j^SSS^iSSSa Similar 
size effect is expected for the Kondo-like scattering on fast TLSs, but we stress that this explanation does not apply 
to the ZBAs in Fig. 122b due to the positive sign of the anomaly. 

B. Zero bias anomalies in metallic glasses 

The zero bias anomalies are also characteristic for metallic glass samples, which were always found to show negative 
anomaly^°'"^°'''^^^^'''^^i^£iiiS22ffi [for examples see Fig.|2] in accordance with the anomalous low-temperature increase 
of the bulk resistanceiSsSiSi^Si In metallic glasses TLSs with a wide distribution of parameters can be present, thus the 
interpretation of the results is even more difficult. Usually the ZBA is attributed to the scattering on TLSs, however 
the magnetic origin also cannot be excluded. The anomalies cannot be suppressed by magnatic field in contrast to the 
magnetic Kondo effect, however the anomalous part of the temperature dependence of the resistance shows a small 
but well defined magnetoresitance, which is not present at higher temperaturesiiisl^£iii££222£ 

Another interesting feature is the telegraph fiuctuation of the ZBA demonstrated by the upper curves in Fig. 
This telegraph noise is only present in the narrow voltage region of the ZBA, that is the amplitude of the fiuctuation 
vanishes as the voltage is increased fFig. l25|l . This extraordinary behavior may be understood by assuming that a slow 
two level fluctuator is not directly changing the resistance (like in Fig. I2UI) . but it is modulating the parameters of a 
fast TLS in the neighborhood. Therefore, the amplitude of the fluctuation only depends on the parameters of a single 
or a few TLSs with the slowly moving defect being in one or the other metastable position, and the contributions 
of electron-phonon interaction, or other TLSs being more far away from this particular defect are less significant. 
Comparing the voltage dependence of the fluctuation with theory, it turns out that the model of fast Kondo-like TLSs 
is in good agreement with the measurement, while the theory of slow TLSs gives a poor fit because of the long tail of 
the PC spectrum.^aSii^i22^ 



C. High frequency behavior of the zero bias anomalies 



The common feature in the yet reported ZBAs is that the anomaly appears in the voltage region ~ 0.1 — 5mV. 
Usually the anomaly shows a negative sign, however even in these cases the elastic contribution from slow TLSs should 
be considered. A direct distinction between the fast and slow TLSs could be made by studying the relaxation dynamics 
of the system. In terms of the two models the position of the singularity reflects the Kondo temperature (Tk) or the 
energy splitting (E), respectively. Inserting the typical voltage range of the anomaly into the models, the standard 



Korringa-like relaxation rate of slow TLSs due to the creation of an electron and hole pair is ~ 10 



10- 



(see 



Eq. I172() , while in the two channel Kondo model the relaxation time is in the range of - io-iis,C°^i999 Therefore the 
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FIG. 24: Differential resistance Rd as a function of bias voltage for FesoB^o (1 and 2) and Fe^2Ni^QCr\4,P\2B(i (3) MCB 
junctions (T = 1.2K). (1) A 6.657 contact displaying almost no noise. (2) A 366f2 contact that shows clear noise around zero 
bias. The noise amplitude decreases as the bias voltage increases. (3) A 145$! contact, showing a clear two-level switching 
behavior between two different Rd peaks (the time dependence of the re sistance at a fix ed bias shows a telegraph noise similar 
to those in Fig. I2UII . Curve 3 has been shifted for clarity. [Taken from lKeiisersl9 96£? with kind permission of the American 
Institute of Physics.] The same paper has also reported about the telegraph fluctuation of the zero bias anomaly in silver break 
junctions. A similar telegraph fluctuation was was observed in the superconducting characteristics of Nia;Nbi_2, metallic glass 
break junctionsi"'^^^""'^''^""" 




FIG. 25: Conductance difference for mechanically controlled metallic glass break junctions. The open squares are obtained 
from the measured resistance change values shown in curve 3 of Fig. 1241 with uncertainties of order 0.1e^//i. The solid line is 
obtained from the theory of Vladar and Zawadowski— — — — - using a Kondo scale of Tk ~ 35K, and two different splittings 
for the fast TLS of E\ = 3.5meV and E2 = 8meV, with the assumption that the modulation by the slow center in duces these 
two different splitting values. The other curve (flUed triangles) represents another measurement from lKeiisersl99'6al . but in this 
case the theoretical flt implies that two fast TLSs are being modulated by a slow one. The inset shows conductance curves 
generated by the theory of Kozub and Kn\\\i^SS^S^ which are also shown in Fig. IIUI in the present review. [Taken from 
IZarandl998l . with kind permission of the American Institute of Physics.] 



boundary between the two cases can be investigated with microwave irradiation measurements in the GHz frequency 
range. 

Experimentally the high frequency measurements were usuaUy performed by irradiating the contact with a chopped 
microwave signal, and measuring the change of the DC current as the microwave is on and off. According to Eq. 11701 
the voltage dependence of the DC shift signal is proportional to the second derivative of the I — V curve if the 
frequency of the irradiation signal is slow compared to the relaxation of the TLS, and the perturbation due to the 
alternating voltage is weak. Previously, it was believed that above the characteristic frequency of the TLS the DC 
shift signal should vanish, however the calculations were only performed for the second generation signal in the case 
of elastic scattering i^Sai^iS^^ Measurements on various systems have shown that the peak in the DC shift signal 
cannot be suppressed in the frequency range of the measurements (< 50 GHz) iSsi^SiSliiiiiSSLEalJsgiiliiSSSi Therefore, 
these studies concluded that the anomaly comes from fast TLSs with a characteristic frequency ujtls ^ 10"'^^ s""'^. 
However, according to our results in Sec. IIV C 41 these measurements are not decisive as the DC shift signal should 
not vanish for slow TLSs either, it should be only slightly reduced, (see upper panels in Fig. El and Fig. [TH|) . The 



49 




bias voltage (mV) 



FIG. 26: RF response signals at 60 GHz (thick solid line) and 0.6 Ghz (dashed line) and d^V/dI^{V) dependence (thin solid 
line) for (a) a 26-Sl Fe8oB2o point contact, (b) a 8-11 Fe78Mo2B2o point contact, and (c) a 15-$! Fe32Ni36Cri4Pi2B6 point 
contact. [Taken from(2steiiiai22& with kind permission of the American Institute of Physics.] 



arbitrariness of the substraction of the background signal makes it even more difficult to draw a clear conclusion 
concerning the frequency dependence of the peak. 

More conclusive results can be obtained by measuring the second harmonic response to the high frequency irradi- 
ation, which is, however, experimentally more difficult. According to our calculations, in the second harmonic signal 
the peak due to the elastic scattering should disappear at high frequency, while the peak due to inelastic scattering 
only decreases by a factor of 2/3 (bottom panels in Fig. 1181 and 1171 respectively). Recent second harmonic gener- 
ation measurements on Ni5gNb4i metallic glasses have shownj^2lJsgiilii222ik that the ZBA cannot be suppressed with 
frequencies up to ~ 5 GHz. This measurement excludes the elastic contribution of slow TLSs due to the lack of 
frequency dependence, and it also excludes the inelastic contribution due to slow TLSs as the sign of the peak is 
negative. However, the authors propose that the frequency dispersion of the smooth background signal might come 
from the scattering on slow TLSs. 

An interesting result has been presented by Balkashin et al. on annealed Ni59Nb4i metallic glass samples i^2iJssiilii222^ 
The DC shift response signal was measured at various microwave irradiation frequencies both before and after an- 
nealing. The unannealed samples have shown insignificant frequency dispersion up the maximum frequency of the 
measurements (^ 5 GHz), in agreement with the previous studies. In the annealed samples, however, the ZBA was 
completely suppressed at high microwave frequencies [Fig. |^. As we have discussed, the ZBA in the DC shift re- 
sponse signal should not dissappear for slow TLSs either, thus the experimental observations on the annealed samples 
cannot be interpreted in terms of the presented theory for slow TLSs. This observation questions whether the ZBA 
originates from the scattering on TLSs in the studied Ni2;Nbi_a; metallic glasses. 

Up to now, the high frequency measurements were only performed on metallic glass samples. Unfortunately in 
metallic glasses various interactions coexist, which makes the interpretation of the results difficult. Especially serious 
problem is the normalization of the measured response signals. In pure samples the curves can be normalized to the 
phonon peaks, as the electron-phonon relaxation time is very small (Toi_ph ~ 10"^'^ s). In metallic glasses, however the 
phonon peaks cannot be identified, and a smooth background is observed, whi ch can contain various contributions. 
In disordered systems the disorder induced electron-electron interaction (see e . g . I Alt shuler 1 99 1|) and the scattering on 
nonequilibrium phonona^^liii^iSSSi^aaSSiiiSS^ are also important, furthermore the scattering on a variety of TLSs with 
different excitation energies can also contribute to the background. Additionally, the materials under study contain 
magnetic ions, thus the interactions of magnetic origin also cannot be excluded. Due to the complexity of the system 
the background signal can also have frequency dispersion in the frequency range of the measurements, which hampers 
the comparison of the response signals measured at different frequencies. 

A more clear picture could be obtained by investigating the high frequency behavior of the ZBAs observed in clean 



50 



1.0 

If) 

£ 0.5 
^. 

-0.5 
-1.0 



1 








M 








1 


\ J 


-30-20 - 10 ( 


) 10 20 30 




meV 



FIG. 27: Evolution of the PC spectra of a contact of the annealed alloy Ni59Nb4i with Ag in measurements at different 
frequencies: a sonic frequency (1) and microwave frequencies ui — 0.48 GHz (2), 2.036 GHz (3), and 4.72 GHz (4). The anomaly 
is strongly suppressed with the microwave irradiation, whereas the ZBA in unannealed samples is unaffected by the same 
irradiation. [Taken from£alkasliin2003 with kind permission of the American Institute of Physics.] 



contacts of pure metals. In this case the ZBA is attributed to the scattering on a single or a few TLSs, and due to 
the baUistic nature of the contact clear phonon peaks are observed at higher energies. In this case the microwave 
response signals are easily compared by normalizing the curves to the phonon peaks, and thus the measurement 
of the second harmonic response signal provides explicit information about the relaxation dynamics of the system. 
Especially interesting would be the study of positive ZBAs in copper nanojunctions (or in other metals like Ag, Au 
and Pt), which are explained by the elastic scattering on slow TLSs. 



D. Conclusions 



A large variety of point contact spectroscopy measurements show ZBAs, which may be attributed to the scattering 
of electrons on TLSs. Certain measurements on Cu nanocontacts show positive ZBAs, which can be well explained by 
the elastic scattering on slow TLS. In this case, however, the reason for the consequently positive sign of the anomaly 
is not well understood. Another group of measurements show negative ZBAs, and among these the experiments 
on Cu nanocinstrictionsSaieiiiSS^ are well interpreted in the framework of the two channel Kondo model. In these 
studies the ZBAs are attributed to a few or even a single TLS, which allows a direct comparison of the experimental 
observations with the theoretical predictions. ZBAs are frequently observed in metallic glass samples, however, in 
this case the complexity of the system makes the interpretation of the results more difficult. Interesting attempts 
were made by investigating the frequency dispersion of the ZBAs in metallic glasses, however some of the results 
were misinterpreted due to the lack of the theoretical calculation for the frequency dependence of the DC shift signal. 
Novel second harmonic generation measurements on Nij;Nbi_3; metallic glasses has proved that in this system the ZBA 
cannot be explained by the scattering on slow TLSs, however slow TLSs can give contribution to the background. 
The study of the high frequency behavior of a few TLSs near clean metallic contacts might provide a more clear 
information about the relaxation dynamics of TLSs. 

In the cases cited above the origins of the TLSs are not well known. TLSs can occur, however, in crystalline material 
with deviation from the stochiometric composition or disorder as well. The TLS is formed in crystalline neighborhood 
where the size or the shape of the cavity in which the atom is sitting can be modified by the surroundings. In 
such cases clear indication of orbital Kondo effect was found in PbTe if some Pb atoms are replaced by smaller Ge 
atoms i^^SiSSiiBiSSi However, this material is not so easy to produce in a controlled way. Recently a new family of 
material was discovered, which are metallic and they are composed of U, Th, As, Se and s,£gaJsi£222i The systematic 
study of these well-controlled materials may also provide the possibility to understand the TLS-electron interaction 
in more detail. 



VI. FINAL REMARKS 



The point contact spectroscopy is a very powerful method to investigate the dynamics of slow TLSs and magnetic 
Kondo impurities and resonances of unknown, nonmagnetic origin, presumably due to some structured defects at zero 
bias. That spectroscopy is complementary of tunneling spectroscopy with oxide barrierJ^"'^^^*^ or break junction in 
the tunneling iiniit,I^SIggSi22SLI^Slg£S222S but there are similarities and differences as well. The local excitations in the 
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barrier as molecular vibrations or magnetic spin flip processes open new tunneling channels resulting in extra current 
at voltages above the excitation energy. That gives just the opposite sign compared to the point contact case, as those 
anomalies are the negative ones. The magnetic impurities in the electrodes nearby the barrier can also depress the local 
tunneling density of states around the Fermi level in the energy range of the Kondo temperatureiMS£ii22ii2SiaSiii2S 
That gives a negative zero bias anomaly just like in the point contact case. 

In the future it would be very important to make systematic studies concerning the TLS, starting with the slow 
TLS with contributions of arbitrary signs due to dominating elastic processes. The detailed studies of the frequency 
dependence based on the theoretical part of the present review would be especially useful and that would clarify some 
of the contradictions of the present time. 

It is very important to study further the resonances of unknown origin which have strong resemblance to nonmag- 
netic Kondo behavioriSasi^SS and to make attempts to determine their origin in at least a few cases. 

The origin of those resonances has been the subject of theoretical debates in the recent years. Originally 
it was proposed that Kondo effect could arise from tunneling TLSs with non-com mutative couplings e.g. due 
to the screening by electrons and electron assisted tunneling (see e.g. ICoxl999l) . It was shown, however, 
that the adiabatic screening of the moving atom by electrons depresses the tunneling rate in such a strong 
manner,^^g'^"iQ^^-^'^'^^"iQ^^-^'""°''^°°i'^i^'""''^^°°i'" that the actual Kondo temperature is negligible. Recently, it has 
been suggested that the situation can be changed if the atom in the double potential well has the ground-state below 
the potential barrier, but the first excited state is above it,^S£^a222^ Then this objection could be avoided, but the 
value of the splitting raises many questions. That theory for an intermediate heavy atom gives too small Tx, not 
more than 0.1 — 0.2 K and that cannot explain the resonances of typical width of 10 — 20 K. There is, however, the 
possibility of light atoms like hydrogen(2ssd2£SJsi2224 as the dominating energy scale is the kinetic energy of the atom 
in the potential well. Therefore, it would be very important in the future to study hydrogen defect either in atomic 
states or water contamination, where the hydrogen could still change its position. 

Summarizing, the origin of slow TLS and the dubious resonances deserve extended further studies. 
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